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Due  to  industrial  development  and  population  growth,  global  demand  for  energy  has  increased 
enormously;  the  increased  consumption  of  primary  sources  of  energy  such  as  coal,  oil  and  natural  gas 
has  exerted  a  strong  influence  on  the  atmospheric  environment.  Among  all  the  alternative  fuels,  hydrogen 
offers  the  greatest  potential  benefits  to  energy  supply  and  the  environment.  However,  combustion  of  pure 
hydrogen  is  challenged  by  the  difficulties  in  production,  storage,  and  end-use.  A  more  reasonable  method 
to  enhance  combustion  is  to  use  hydrogen  as  an  additive  to  fossil  fuels. 

This  paper  presents  a  comprehensive  review  of  recent  progress  in  fundamental  investigations  of 
hydrogen  enriched  combustion.  A  broad  range  of  topics  on  fundamental  combustion  of  hydrogen/ 
hydrocarbon  mixtures  are  discussed.  The  minimum  ignition  energy  of  hydrogen/hydrocarbon  blends 
shows  significant  decrease  with  the  increased  hydrogen  fraction  XH2.  With  an  increase  of  the  hydrogen 
fraction,  ignition  behavior  of  hydrogen/hydrocarbon  fuel  blends  exhibits  three  distinct  behaviors. 
Additionally,  the  laminar  flame  speed  of  hydrogen/hydrocarbon  fuel  blends  is  found  to  exhibit  three 
behavior  regimes.  Furthermore,  the  flame  front  instability  behavior  is  influenced  by  the  hydrogen  addition 
through  the  combined  effect  of  body  force,  the  hydrodynamic  effect  and  the  thermal-diffusion  effect. 
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1.  Introduction 

Fossil  fuels  in  the  form  of  coal,  oil,  and  natural  gas  have  been 
used  in  transportation,  heating,  and  power  plants  since  the 
beginning  of  the  industrial  revolution  in  the  18th  century.  How¬ 
ever,  increasing  global  demand  for  energy  due  to  industrial 
development  and  population  growth  is  being  threatened  by  the 
limited  supply  of  fossil  fuel.  Additionally,  the  fact  that  oil  is  found 
primarily  in  only  a  few  regions  of  the  world  creates  a  strong 
economic  dependence,  especially  in  developing  countries.  Further¬ 
more,  the  combustion  of  fossil  fuel  produces  carbon  dioxide, 
which  has  been  named  as  a  primary  cause  of  global  warming; 
and  the  chemical  and  particulate  emissions  produced  from  burn¬ 
ing  of  fossil  fuel  put  human  health  at  risk. 

3.2.  The  status  of  hydrogen  technology:  production,  storage, 
and  end-use 

Research  on  alternative  fuels  have  been  driven  by  environ¬ 
mental  and  energy  supply  challenges.  Among  the  various  alter¬ 
native  fuels,  hydrogen  offers  the  significant  potential  benefits  to 
the  environment  and  energy  supply  [1],  Hydrogen  is  the  most 
abundant  element  on  Earth,  but  only  1%  exists  in  the  form  of 
molecular  hydrogen  H2.  There  are  a  variety  of  widely  available 
methods  for  producing  hydrogen.  Currently,  commercially  avail¬ 
able  hydrogen  is  primarily  produced  through  steam  reforming 
from  natural  gas,  in  which  carbon-fuel  is  still  used  [2], 

Cost  effective  ways  to  mass  produce  hydrogen  from  non-carbon 
fuel  is  a  challenge  to  the  use  of  hydrogen  as  a  substitute  for  fossil 
fuels.  The  most  promising  method  is  splitting  water.  However,  it 
takes  energy  to  split  water  molecule,  even  though  that  energy 
is  later  replaced  by  energy  released  from  the  oxidation  of 
hydrogen  [3],  Non-carbon  sources  of  energy  such  as  solar,  nuclear, 
hydroelectric,  or  wind  energies  have  been  tested  in  the  water¬ 
splitting  process,  the  goal  being  to  eliminate  fossil  fuels  from  the 
water-splitting  cycle. 

Regarding  the  issue  of  on-board  storage,  a  US  Department  of 
Energy  (DOE),  Hydrogen  Program  has  set  a  standard  [4]:  the 
storage  system  should  contain  a  hydrogen  weight  of  at  least  6.5% 
of  the  system  weight,  and  the  volumetric  density  of  62  kg  H2/m3. 
Conventionally,  hydrogen  is  stored  in  a  high  pressure  steel  tank, 
with  a  typical  filled  pressure  of  up  to  200  bar  in  most  countries. 
Carbon/polymer  tanks,  or  fiberglass/aluminum  tanks,  can  with¬ 
stand  pressures  as  high  as  600  bar,  but  the  volumetric  energy 
density  is  still  low.  Storage  of  H2  on  activated  carbon  (AC) 
materials  at  temperatures  below  150 1<  has  been  investigated 
and  proved  unsatisfactory  [5].  Liquefied  hydrogen  has  better 
volumetric  energy  density,  but  the  process  of  liquefaction  requires 
large  amounts  of  energy  to  create  an  environment  with 


temperatures  as  low  as  20  1<  and  perfect  insulation  techniques. 
Some  metals  or  alloys  can  form  hydrides,  which  can  release 
hydrogen,  but  the  gravimetric  density  is  typically  less  than  3% 
due  to  the  heaviness  of  the  metals.  Single-walled  carbon  nano¬ 
tubes  (SWNTs)  can  absorb  large  amounts  of  hydrogen  into  its 
molecularly  dimensioned  pores.  Fig.  1  shows  that  none  of  the 
hydrogen  storage  techniques  is  capable  of  meeting  DOE  goals  [6], 
but  there  is  reason  to  expect  that  one  day  better  hydrogen  storage 
materials  and  techniques  will  be  discovered  and  developed. 

There  are  essentially  two  methods  for  hydrogen  usage:  hydro¬ 
gen  fuel  cells  and  hydrogen  internal  combustion  engines.  A  fuel 
cell  converts  the  chemical  energy  of  the  fuel  directly  to  electricity. 
A  hydrogen-oxygen  fuel  cell  has  an  ideal  electrical  conversion 
efficiency  of  83%.  In  a  real  fuel  cell,  up  to  60%  of  the  chemical 
energy  is  converted  to  electricity,  and  the  remainder  is  converted 
to  heat.  Theoretically  many  other  fuels  can  be  used  in  fuel  cells, 
but  hydrogen  is  the  best  candidate  because  of  the  rapid  reaction 
kinetics  at  the  anode.  Worldwide  efforts  have  been  made  to 
commercialize  the  proton  exchange  membrane  (PEM),  fuel  cell 
in  the  last  few  years;  however,  the  cost  of  PEM  fuel  cells  is  still 
much  higher  than  conventional  internal  combustion  engines. 
The  cost  of  power  for  PEMs  fuel  cells  for  transportation  is 
$1500-10,000  per  kilowatt,  while  the  power  cost  of  conventional 
automobile  is  less  than  $50  per  kilowatt  [3], 

1.2.  Hydrogen  as  an  enriching  agent 

Lean  combustion  of  hydrocarbon  in  nearly  all  combustors  such 
as  gas  turbines,  furnaces,  boilers  and  internal  combustion  engines 


Fig.  1.  Installed  energy  densities  for  several  vehicular  hydrogen  storage  techniques, 
from  Ref.  [6], 


C.  Tang  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  30  (2014)  195-2 16 


197 


Nomenclature 

Ead 

adiabatic  flame  temperature 

V 

volume 

4> 

equivalence  ratio 

le 

Zeldovich  number 

a 

thermo  expansion  rati 

K 

stretch  rate 

Abbreviations 

Ea 

activation  energy 

El,  min 

minimum  ignition  energy 

CCMHDI 

combustion  chemistry  of  methane  and  hydrogen 

I(a 

Karlovitz  number 

dominated  ignition 

P 

pressure 

HCDI 

hydrogen  chemistry  dominated  ignition 

Rh 

hydrogen  addition  parameter,  see  Eq.  (11) 

LPC 

lean  premixed  combustion 

T 

temperature 

MIE 

minimum  ignition  energy 

U 

internal  energy 

PR 

pressure  rise 

•Xh2 

hydrogen  fraction  in  the  fuel  blends 

UFL 

upper  flammability  limit 

4'r 

Fuel  based  equivalence  ratio,  see  Eq.  (10) 

HCNG 

hydrogen/compressed  natural  gas 

Kext 

extinction  strain  rate 

LFL 

lower  flammability  limit 

dq 

quenching  distance 

MCDI 

methane  chemistry  dominated  ignition 

k 

linearity  coefficient,  see  Eq.  (13) 

NG 

natural  gas 

f 

mass  burning  flux 

SFHC 

stoichiometric  fraction  of  the  hydrocarbon  in  the  fuel 

Le 

Lewis  number 

mixtures 

n 

flame  radius 

VI 

visualization 

5° 

laminar  flame  speed  Tu  unburned  gas  temperature 

is  advantageous  in  terms  of  NO*,  CO  and  unburned  HC  emissions 
as  well  as  for  fuel  economy.  However,  lean  hydrocarbon/air 
mixtures  are  hard  to  ignite  and  have  low  flame  speed,  which  can 
potentially  cause  misfire  and  results  in  increased  unburned  hydro¬ 
carbons,  reduced  performance  and  poor  fuel  economy.  In  addition, 
lean  hydrocarbon  flame  is  hard  to  burn  steadily,  which  causes 
malfunction  of  the  combustors. 

The  properties  that  contribute  to  use  of  hydrogen  in  internal 
combustion  engines  are  its  wide  flammability  range,  which 
enables  hydrogen  to  run  on  a  very  lean  mixture;  low  ignition 
energy,  which  ensures  prompt  ignition;  high  reactivity  and  flame 
speed,  which  moves  the  real  engine  cycle  closer  to  the  thermo¬ 
dynamic  ideal;  high  diffusivity,  which  facilitates  the  uniformity  of 
fuel  and  air;  and  high  auto-ignition  temperature,  which  allows 
larger  compression  ratio  of  the  engine.  However,  when  pure 
hydrogen  is  used,  the  engines  are  prone  to  pre-ignition  due  to 
hydrogen’s  lower  ignition  energy,  wider  flammability  range  and 
shorter  quenching  distance.  In  addition,  when  pure  hydrogen  is 
burned  in  the  engine,  the  heat  release  is  very  fast,  and  the 
temperature  is  high;  thus  the  engine  could  potentially  emit  high 
levels  of  NOx  and  engine  reliability  is  challenged  because  of  the 
high  thermal  load. 

To  some  extent  the  drawbacks  of  burning  pure  hydrogen  or 
pure  hydrocarbons  can  be  reduced  simultaneously  by  mixing 
hydrocarbon  with  a  certain  amount  of  hydrogen.  The  low  ignition 
energy  of  hydrogen  makes  hydrogen/hydrocarbon  mixtures  easier 
to  ignite.  Hydrogen’s  high  flame  speed  reduces  misfire,  thereby 
reducing  emissions.  The  energy  density  of  the  mixture  is  augmen¬ 
ted  because  of  the  increased  hydrogen-to-carbon  ratio;  thus  the 
torque  is  improved  at  wide  open  throttle  conditions.  Using 
hydrogen  as  a  hydrocarbon  combustion  stimulant  has  received 
considerable  attention  as  a  near  term  solution,  before  a  full 
hydrogen  infrastructure  -  in  terms  of  production,  storage  and 
end-use  -  is  established. 

1.3.  Outline  of  the  review 

Following  is  a  review  of  the  wide  spectrum  of  research 
activities  on  hydrogen  stimulated  hydrocarbon  combustion.  I  In 
Section  2,  the  flammability  limits  of  hydrogen  enriched  hydro¬ 
carbons  are  discussed.  In  Section  3,  the  literature  on  the  effect  of 


hydrogen  addition  on  spark  ignition  characteristics  is  presented, 
followed  by  a  discussion  on  the  hydrogen  addition  effect  on  auto¬ 
ignition  by  shock  waves,  in  Section  4.  The  hydrogen  addition 
effects  on  laminar  flame  speed  and  flame  front  instability  of 
hydrocarbon/air  mixture  are  reviewed  in  Section  5.  In  Section  6, 
the  hydrogen  addition  effect  on  lean  premixed  combustion  stabi¬ 
lity  and  emissions  from  lean  premixed  combustion  is  reviewed. 

2.  Flammability  limits 

A  lower  flammability  limit  ( LFL )  and  an  upper  flammability 
limit  (UFL)  are  typically  determined  for  a  fuel-oxidizer  system  at  a 
given  pressure  and  temperature,  and  when  the  fuel/oxidizer  ratio 
is  lower  than  LFL,  or  higher  than  UFL,  the  mixture  is  inflammable. 
The  studies  on  flammability  issues  are  primarily  motivated  by  its 
importance  in  assessing  fire  hazards.  In  addition,  advantages  in 
lean  combustion,  such  as  improved  emissions  and  efficiency 
characteristics,  are  limited  by  the  LFL. 

2.1.  Flammability  limits  of  hydrocarbons 

There  are  several  experimental  techniques  which  determine 
flammability  limits.  The  most  direct  and  best  known  method  is  the 
so-called  visual  identification  method  (VI),  developed  by  the  former 
US  Bureau  of  Mines.  In  this  method,  the  mixture  is  defined  as 
“flammable"  if  the  flame,  induced  by  a  spark  at  the  bottom  of  a 
vertically  oriented  glass  tube,  can  propagate  upwardly  a  certain 
distance.  Much  of  the  work  with  this  method  was  performed  and 
summarized  by  Coward  and  Jones  [7],  Zabetakis  [8],  and  Kuchta  [9], 
Table  1  shows  some  examples  of  the  flammability  limits,  determined 
using  the  visual  identification  (VI)  method.  It  is  seen  that  hydrogen 
has  the  lowest  LFL,  indicating  that  hydrogen  is  potentially  capable  of 
burning  under  very  lean  mixture  conditions. 

Another  apparatus  for  measuring  the  flammability  limit  is  a 
closed,  spherical  reaction  vessel  with  central  ignition.  In  this  method, 
the  flammability  criterion  is  the  relative  pressure  increase  from 
combustion.  The  Pittsburgh  Research  Laboratory  has  published 
considerable  flammability  data  using  the  pressure  rise  (PR)  criterion 
[10-12],  A  counter-flow  burner  can  also  be  used  to  determine  the 
flammability  limit.  In  this  method,  the  global  strain  rate  (a  gradient  of 
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Table  1 

Flammability  limits  of  some  common  fuels  at  1  atm,  25  C.  Data  from  Ref.  [8], 


Fuel 

<p 

LFL 

UFL 

Hydrogen 

0.1 

7.14 

Carbon  monoxide 

0.34 

6.8 

Ammonia 

0.63 

1.4 

Methane 

0.5 

1.67 

Ethane 

0.52 

2.4 

Propane 

0.56 

2.7 

Butane 

0.57 

2.8 

Ethylene 

0.40 

8.0 

Acetylene 

0.31 

oo 

Benzene 

0.56 

3.7 

Methanol 

0.51 

4.0 

Ethanol 

0.41 

2.8 

the  exit  velocity)  is  plotted  as  a  function  of  fuel  concentration.  The 
flammability  limit  can  be  obtained  by  extrapolating  the  curve  to  the 
minimum  strain  rate  [13]. 

Various  organizations  have  attempted  to  standardize  the  mea¬ 
surements  of  flammability  limits.  To  list  some  examples:  The 
American  Society  for  Testing  and  Materials  (ASTM)  has  used  a 
5L  glass  sphere  (ASTM  E681),  to  visually  detect  the  flame 
propagation  in  order  to  define  the  flammability  [14],  a  4L  near- 
spherical  vessel  (ASTM  E2079),  with  high  ignition  energy  and  the 
7%  pressure  rise  criterion  [15],  and  a  1  L  minimum  vessel  with 
76  mm  diameter  (ASTM  E918),  with  7%  pressure  rise  criterion, 
capable  of  operating  under  elevated  temperatures  and  pressures 
[16].  The  European  Standard  EN  1839  is  subdivided  into  EN  1839T 
and  EN  1839B.  EN  1839T  defines  flammability  by  detecting  10  cm 
of  flame  propagation  upon  ignition  in  a  glass  cylinder,  and  the  EN 
1839B  uses  the  5%  pressure  increase  in  a  cylinder  or  spherical 
vessel  with  volume  larger  than  5  L.  ISO  10156  recommended  ASTM 
E181  and  EN  1839T  as  the  experimental  criteria  for  flammability 
determination.  Van  den  Schoor  et  al.  [17]  and  Schroder  and 
Daubitz  [18]  compared  EN  and  ASTM  methods  and  their  investiga¬ 
tion  showed  similar  results  for  the  most  important  gases  and  gas 
mixtures. 

Increasing  computational  capability  and  development  of  reac¬ 
tion  mechanisms  makes  possible  the  numerical  simulation  of 
limited  flames  with  different  geometries.  Simulation  attempts 
include  ID  planar  flame  [19,20]  and  the  quasi-ID  spherically 
propagating  flame  [21],  with  or  without  radiation  heat  loss. 
Theoretically,  there  exists  a  limiting  flame  speed  below  which 
the  flame  is  unable  to  propagate;  this  approach  had  been  used  to 
determine  flammability  limits  [22,23].  Analogous  to  the  limiting 
flame  speed  approach,  Burgess  and  Wheeler  [24]  proposed  a 
temperature  limiting  concept,  which  states  that  for  many  lean 
limit  hydrocarbons-air  mixtures,  the  LFL  (in  terms  of  fuel  volume 
fraction)  times  the  heating  value  of  the  fuel,  <jci,  is  nearly  a 
constant  around  920  kcal/mole,  limiting  the  temperature  to 
around  1450  K  [25,26].  This  concept  is  implicitly  used  when 
applying  Le  Chatelier’s  mixing  rule,  which  is  based  on  the 
assumption  that  the  heat  of  combustion  of  the  mixture  at  the 
lean  flammability  limit  is  the  heat  of  combustion  of  its  composing 
fuels  at  their  respective  LFLs,  as  shown  in  the  following  equation: 


(LFL)mix  = 


X, 


i  =  1  (LFL),- 


(1) 


where  X,  is  the  volume  fraction  of  the  fuel  i  and  (LFL),  is  its 
corresponding  flammability  limit,  and  n  is  the  total  number  of  fuel 
components. 


2.2.  Flammability  limit  of  hydrogen/hydrocarbon  mixtures 

Several  researchers  have  tested  the  effects  of  hydrogen  addition 
on  the  flammability  limits  of  premixed  hydrocarbon  air  mixtures 
with  a  variety  of  methods  described  in  Section  2.1.  Wierzba  and 
Ale  [27]  investigated  the  effect  of  hydrogen  addition  on  the  rich 
flammability  limit  of  methane,  ethylene,  and  propane  at  elevated 
temperatures  with  the  VI  method.  Their  measured  UFL  signifi¬ 
cantly  increases  when  hydrogen  is  added,  and  for  methane  and 
propane  the  measurements  agree  well,  but  for  ethylene  the  results 
deviate  remarkably  with  the  predictions  by  Le  Chatelier’s  Rule. 
Shoshin  and  de  Goey  [28]  also  measured  the  flammability  limit  of 
methane/hydrogen/air  mixtures  with  the  VI  method  -  especially 
the  effects  of  tube  diameter.  Miao  et  al.  [29]  measured  the 
flammability  limit  of  natural  gas/H2  and  methane/H2  blends  using 
the  PR  method;  their  results  showed  no  difference  between 
natural  gas  and  methane.  In  Ref.  [17],  Van  den  Schoor  et  al.  applied 
the  VI  and  PR  method,  as  well  as  a  planar  and  spherical  flame 
model  to  a  methane/hydrogen/air  system;  the  limiting  burning 
velocity  of  5  cm/s  criterion  predicts  the  UFL  well,  as  compared 
with  the  experimental  results.  Guo  et  al.  [30]  studied  the  hydrogen 
addition  effect  on  the  flammability  limit  of  an  ultra-lean  methane/ 
air  stream.  As  shown  in  Fig.  2,  the  upper  branch  is  the  stretch 
extinction  limit,  caused  by  high  stretch  rate,  and  the  lower  branch 
is  the  radiation  extinction  limit,  due  to  an  ultra  low  stretch  rate.  It 
is  seen  that  with  the  increase  of  XH2,  both  limits  shift  to  a  leaner 
mixture,  thus  the  turning  point  of  the  “C”  shaped  curve,  defined  as 
the  LFL,  is  consequently  decreased. 

The  recent  literature  reporting  flammability  limit  data  of 
hydrogen/methane/air  mixtures  was  summarized  in  a  flammabil¬ 
ity  regime  diagram  in  the  parameter  space  of  hydrogen  fraction 
XH2  and  equivalence  ratio  <p,  as  shown  in  Fig.  3.  The  flammable 
and  nonflammable  regimes  are  separated  by  the  LFL  and  UFL. 
The  experimentally  measured  LFLs  with  different  techniques  are 
consistent.  For  the  UFLs,  large  discrepancies  exist  among  different 
literatures.  The  predictions  by  Le  Chatelier’s  Rule  were  also  shown 
in  Fig.  3;  data  of  LFL  and  UFL  of  methane  and  hydrogen  in  Table  1 
were  used.  For  LFL,  Le  Chatelier’s  prediction  is  in  good  agreement 
with  those  measurements,  but  for  UFL,  the  prediction  deviated 
considerably  from  most  of  the  experimental  measured  results. 


3.  Hydrogen  assisted  ignition 

The  flame  is  initiated  through  two  methods:  one  is  forced  ignition 
by  an  electric  spark,  laser,  hot  air  jets  or  plasma;  the  other  is  auto¬ 
ignition  by  shock  waves,  which  will  be  discussed  in  Section  4. 


* 


Fig.  2.  Extinction  strain  rate  as  a  function  of  equivalence  ratio  for  counter-flowing 
methane/hydrogen/air  mixtures  [30]. 
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Fig.  3.  Flammability  region  diagram  for  methane/hydrogen/air  mixture  at  1  atm 
and  room  temperature  [17,27-30]. 


Fig.  4.  Correlation  of  the  MlEs  with  quenching  diameter  for  different  fuels.  Adapted 
from  Calcote  et  al.  [33], 

3.1.  Spark  ignition  of  premixed  hydrogen/hydrocarbon/air  mixtures 

3.1.1.  Thermal-diffusion  theory  on  M1E 

Generally,  ignition  of  combustible  mixture  is  initiated  by  an 
electric  spark,  and  the  subsequent  flame  is  propagated  outward  in 
an  expanding  sphere.  The  literatures  show  that  the  spark  ignition 
progresses  in  three  stages:  the  breakdown  stage,  the  arc  stage  and 
the  glow  stage.  The  breakdown  stage  occurs  in  a  timescale  of 
10~6s,  during  which  a  conducting  channel,  consisting  of  high 
temperature  equilibrium  plasma,  is  opened.  In  the  subsequent  arc 
stage,  the  plasma  kernel  grows  due  to  the  deposit  of  electric 
energy.  The  last  stage  is  the  ignition  of  the  combustion  mixture 
surrounding  the  hot  plasma  kernel.  The  overall  energy  deposit  by 
the  spark  is  the  most  important  factor  for  successful  ignition. 
Minimum  ignition  energy  (MIE)  is  the  minimum  energy  required 
to  ignite  a  combustible  mixture.  In  the  thermal-diffusion  theory 
[31,32],  MIE  has  been  constantly  correlated  to  the  quenching 
distance,  dq,  or  laminar  flame  thickness,  <5|.  For  fuels  with  different 
molecular  structures,  Calcote  et  al.  [33]  have  correlated  their  MIEs, 
£i,min.  with  the  measured  dq  as  shown  in  Fig.  4.  The  slope  of  the 


line  in  the  logarithmic  scale  yields  the  following  relation: 

E;,min~c£48  (2) 

Williams  [32]  also  empirically  correlated,  and  a  power  of  2  was 
proposed.  Classman  [34]  and  Law  [25]  stated  that  MIE  should  be 
proportional  to  the  amount  of  energy  needed  to  heat  a  spherical 
volume  of  the  mixture  of  radius  of  the  order  of  laminar  flame 
thickness,  St  (of  the  order  of  quenching  distance),  to  the  burned 
temperature  (close  to  Tad),  yielding 

El,  min  ~  Pu^l  ad —  Tu)  ~  dq  (3) 

To  the  authors’  knowledge,  the  only  literature  which  reports  on 
the  minimum  ignition  energy  and  quenching  distance  for  hydrogen/ 
methane/air  mixture  is  by  Fiankinson  et  al.  [35]  and  Fukuda  et  al.  [36], 
The  minimum  ignition  energy  and  the  quenching  distance  both 
decrease  with  the  increase  of  XH2,  indicating  that  the  addition  of 
hydrogen  is  potentially  capable  of  facilitating  successful  ignition. 

Fi.min  and  dq  for  methane/hydrogen/air  mixtures  were  corre¬ 
lated  in  Fig.  5  for  hydrogen  fractions  of  0,  25%,  50%,  75%  and  100%. 
The  data  were  obtained  directly  or  through  interpolation  from 
Refs.  [35,36],  It  is  seen  that  Eimin  and  dq  for  all  mixture  conditions 
could  be  linearly  correlated  in  the  logarithmic  scale,  yielding 

E,.  min  ~d2q32  (4) 

in  which  the  power  of  2.32  compares  favorably  with  that  deter¬ 
mined  by  Calcote  et  al.  [33],  Eqs.  (2)  and  (4)  are  obtained  by  fitting 
the  experimental  data,  while  Eq.  (3)  is  derived  based  on  the 


cfq  /  mm 

Fig.  5.  Minimum  ignition  energy  as  a  function  of  quenching  distance  for  methane/ 
hydrogen/air  mixture  at  1  atm  and  room  temperature.  Data  obtained  from  Refs.  [35,36], 


Fig.  6.  Rate  of  production  of  H20  profile  for  XH2  <  20%.  Adapted  from  Ju  and  Niioka 

[45]. 
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thermal-diffusion  theory  and  the  level  of  agreement  between  the 
experimental  measurement  and  theoretical  prediction  is  acceptable. 


3.1.2.  Self-sustained  flame  propagation 

The  minimum  ignition  energy  has  been  extensively  measured 
[31,33],  based  on  which  ignition  process  is  better  realized;  how¬ 
ever,  the  flame  dynamics  at  the  initial  stage  were  not  included  in 
these  studies.  Because  the  flame  kernel  radius  was  very  small,  the 
stretch  effect  could  be  large,  and  could  possibly  result  in  failure  to 
form  flame  propagation.  Law  [25]  proposed  a  generalized  factor  to 
quantify  the  transient  effect  of  stretch 

a=Z^(Le-'-\)Ka  (5) 

where 


Ze  = 


Ea(Tad  —  Tu) 
R°TL  ’ 


Ka  =  -wK 


and  k  = 


2  drf 
Tf  dt 


are  respectively  the  Zeldovich  number,  Karlovitz  number  and  the 
instantaneous  stretch  for  outwardly  propagating  flames.  Furthermore, 
Le,  Ea,  Tad,  Tu,  R°,  su°  and  rf  are  respectively  the  Lewis  number, 
activation  energy,  adiabatic  temperature,  unburned  gas  temperature, 
universal  gas  constant,  laminar  flame  speed  and  flame  radius. 

it  is  noted  that  the  Karlovitz  number,  which  demonstrates  the 
effect  of  stretch,  and  the  Lewis  number,  which  represents  the 
effect  of  non-equidiffusion,  are  intrinsically  coupled,  and  le 
represents  the  chemical  weighing  in  Eq.  (5).  In  addition,  outwardly 
propagating  spherical  flame  experiences  positive  stretch.  Eq.  (5) 
shows  that  the  temperature  (hence  the  flame  propagation)  is 
promoted  if  Le  <  1,  while  the  opposite  holds  for  Le  >  1.  In  the 
latter  case,  the  competition  between  the  electric  energy  supply 
and  the  stretch  effect  may  result  in  failure  of  a  self-sustained 
propagating  flame  if  the  flame  radius  does  not  grow  large  enough 
beyond  which  the  stretch  effect  could  be  negligible. 

Recently,  He  [37]  and  Chen  and  Ju  [38],  numerically  investi¬ 
gated  the  propagation  of  spherical  flames,  allowing  for  continuous 
supply  of  ignition  energy.  The  results  showed  that  for  Le  <  1 
mixture,  the  initial  high  flame  speed  due  to  spark  energy  supply 
decreases  to  planar  value  as  the  flame  expands  (provided  that  a 
flame  kernel  is  successfully  established).  While  for  Le  >  1  mixture, 
the  response  is  quite  different.  If  the  ignition  energy  is  sufficiently 
large,  the  flame  speed  decreases  to  a  minimum  value  and  then 
increases  to  the  planar  value.  Otherwise,  the  continuously  decel¬ 
erating  flame  will  extinguish  at  a  certain  flame  radius.  Kelley  et  al. 
[39]  then  experimentally  examined  the  above  theoretical  predic¬ 
tions  with  hydrogen/air  and  butane/air  flames.  Indeed,  the  flame 
propagation  behavior  for  lean  hydrogen/air  mixture  and  rich 
butane/air  mixture  (Le<l)  is  quite  similar,  so  that  the  flame 
propagation  speed  decreases  to  an  asymptotic  value.  However, 
for  the  case  of  very  rich  hydrogen/air  or  lean  butane/air  mixture 
(Le  >  1 ),  the  flame  speed  first  decreases  to  a  minimum,  and  the 
self-sustained  flame  propagation  is  subsequently  recovered. 

Although  there  have  been  no  publications  about  flame  kernel 
growth  for  a  hydrogen/hydrocarbon/air  mixture,  it  is  anticipated  that 
adding  hydrogen  to  a  lean  hydrocarbon/air  mixture  will  influence  the 
flame  kernel  evolution,  because  hydrogen  will  affect  the  preferential 
diffusion  characteristics  of  hydrocarbons-air  mixture,  as  will  be 
discussed  in  detail  in  Section  5.2. 


3.2.  Ignition  of  methane/hydrogen  and  air  in  supersonic  mixing  layer 

The  practical  objective  of  the  studies  on  ignition  in  supersonic 
mixing  layers  is  the  determination  of  ignition  distance  in  response  to 
the  interest  in  the  development  of  scramjet  engines.  The  ignition 
kinetics  of  hydrogen/air  [40,41]  and  methane/air  [42]  has  been 
studied  and  methane  was  found  to  have  a  long  ignition  delay  time. 


Thus  ignition  improvement  is  required.  Silane  is  an  effective  candi¬ 
date,  but  its  highly  reactive  nature  renders  on-board  use  difficult  and 
dangerous  [43],  Hence  hydrogen  could  be  more  practically  used  as 
the  ignition  stimulus.  For  instance,  Bier  et  al.  [44]  conducted 
experiments  on  the  ignition  of  transversely  injected  methane  and 
hydrogen  into  parallel  supersonic  air  stream  and  found  that  the 
lowest  static  air  temperature  for  stable  combustion  of  hydrogen  was 
700  C  -  much  lower  than  1300  C  -  which  is  the  corresponding 
value  for  methane.  However,  there  have  been  only  a  few  investiga¬ 
tions  on  the  ignition  of  methane/hydrogen  vs.  air  stream  in  super¬ 
sonic  mixing  layer.  Ju  and  Niioka  [45]  numerically  investigated  the 
effect  of  hydrogen  fraction  in  the  methane/hydrogen  fuel  stream 
(here,  also  designated  as  XH2)  on  the  ignition  process.  Velocities  of 
the  air  and  fuel  mixture  streams  were  set  to  be  1200  and  1000  m/s, 
respectively.  The  computational  domain  is  16  x  3  cm2.  They  found 
that  the  ignition  behavior  strongly  depends  on  XH2.  As  shown  in 
Fig.  6,  for  XH2  less  than  9%  (although  strong  ignition  was  not  observed 
in  the  computational  domain),  the  concentration  of  the  produced 
H20  still  increased  with  increasing  XH2-  When  XH2  was  further 
increased  up  to  12%,  a  sharp  increase  in  the  rate  of  production  of 
H20  was  observed  in  the  computational  domain,  which  is  the 
manifestation  of  strong  ignition. 

For  a  very  high  hydrogen  fraction  up  to  85%,  the  mixture  behaves 
analogously  to  hydrogen  and  the  presence  of  methane  acts  as  an 
ignition  inhibitor.  Ju  and  Niioka  [45]  identified  three  sub-regimes  with 
increasing  methane  fraction  XCH4  (decreasing  XH2).  As  shown  in  Fig.  7, 
in  the  first  sub-regime,  the  presence  of  even  a  very  small  amount  of 
methane  significantly  increased  the  ignition  distance.  This  significant 
increase  leveled  off  in  the  second  sub-regime.  In  the  third  sub-regime 
as  methane  concentration  was  further  increased,  a  second  rapid 
increase  in  ignition  time  was  observed.  Specifically,  for  XCH4  <  2%, 
the  ignition  time  increased  quickly  with  increasing  XCH4.  In  this  sub¬ 
regime,  the  endothermic  reaction  CH4<->CH3+H  and  the  H  radical 
scavenging  reaction  H+CH4<->CH3+H2  were  both  promoted  with 
more  CH4  in  the  mixture,  resulting  in  a  dramatic  slowdown  of  the 
ignition.  In  the  second  sub-regime,  new  H  radical,  produced  through 
CH4<  >CHj+H,  became  effective,  which  leveled  off  the  rapid  increase 
in  ignition  time.  As  Xch4  was  further  increased,  the  increased  H  radical 
induced  a  fast  procession  of  H + CH4  <->  CH3 + H2,  which  suppressed  the 
main  chain-branching  reaction  H+02<->0H+0,  leading  to  a  second 
rapid  increase  in  ignition  time. 

3.3.  Ignition  temperature  of  hydrogen/methane  by  counter-flowing 
hot  air  jets 

Non-premixed  ignition  in  counter-flowing  fuel,  such  as  H2 
[46-51],  methane  [52],  and  higher  hydrocarbons  [53],  vs.  hot  air 


Fig.  7.  Ignition  distance  as  a  function  of  XH2  for  XH2  >  85%.  Adapted  from  Ju  and 
Niioka  [45],  Three  sub-regimes  were  identified:  (a)  chain-branching  inhibition 
regime,  (b)  transition  regime,  and  (c)  reaction  competition  regime. 
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jets,  has  been  extensively  investigated.  However,  there  has  been 
only  one  paper  on  the  subject  of  ignition  characteristics  of 
hydrogen/hydrocarbon  fuel  mixtures:  Fotache  et  al.  [54]  measured 
the  ignition  temperature  of  counter-flowing  H2/N2  vs.  hot  air  at 
1  atm,  and  the  strain  rate  was  fixed  at  150  /s.  N2  was  then 
substituted  with  CH4  and  the  corresponding  ignition  temperature 
was  also  measured.  Numerically,  the  ignition  temperature  of  H2/N2 
was  calculated  using  the  mechanism  of  Kim  et  al.  [55],  The  GR1 
Mech  1.2  [56]  was  adopted  for  the  H2/CH4  system.  As  shown  in 
Fig.  8,  the  measured  and  the  calculated  ignition  temperature  for 
both  H2/N2  (the  squares  and  the  dotted  lines)  and  H2/CH4  system 
(the  circles  and  the  solid  lines)  agree  well,  except  for  XH2  >  30%  in 
the  H2/CH4  system;  the  GRI  Mechl.2  results  deviated  substantially 
with  the  measurements.  A  modified  mechanism  based  on  GRI 
Mechl.2  was  developed  (dashed  lines).  Three  distinct  regimes 
were  observed.  When  XH2  is  less  than  6-7%,  the  H2/N2  system 
failed  to  ignite  and  the  ignition  temperature  of  H2/CH4  was 
significantly  decreased,  even  with  small  hydrogen  concentrations. 
ForXH2  between  7%  and  30%,  the  ignition  temperature  of  methane 
still  decreased,  but  much  less  significantly  than  in  the  first  regime. 
When  XH2  was  larger  than  30%,  the  ignition  temperature  was 
almost  independent  of  the  hydrogen  concentration. 

In  Section  3.2,  the  work  of  Ju  and  Niioka  [45]  is  reviewed;  it  also 
identified  three  ignition  regimes,  but  at  high  hydrogen  concentra¬ 
tions.  Ju  and  Niioka’s  results  show  the  strong  inhibitive  effect  of 
methane  at  high  hydrogen  concentrations.  This  is  in  contrast  to 
the  results  of  Fotache  et  al.  [54],  in  which  as  much  as  70%  methane 
concentration  changes  ignition  temperature  very  little.  This  qua¬ 
litatively  different  behavior  was  explained  by  Fotache  et  al.  as 
follows:  first,  the  characteristic  residence  time  in  Ref.  [45]  is 
considerably  shorter  than  in  Ref.  [54],  Fotache  et  al.  found  that 
at  high  strain  rates  (104-105  s  1 ),  the  ignition  temperature 
becomes  sensitive  to  methane  addition  -  even  for  a  small  amount 
of  methane  [54],  Second,  the  mechanism  Ju  and  Niioka  used 
imposed  a  higher  rate  of  H  scavenging  reaction,  which  renders  a 
higher  sensitivity  of  methane  addition.  Finally,  the  short  residence 
time  in  Ref.  [45]  requires  the  heating  of  the  fuel  jet  to  1000  K;  thus 
the  reaction  zone  is  closer  to  the  fuel  side  than  in  Ref.  [54],  and  the 
effect  of  methane  is  further  enhanced. 


4.  Shock  tube  ignition  delay  times  of  hydrogen/hydrocarbons 

It  is  widely  recognized  that  knock  in  ICEs  is  caused  by  the  auto¬ 
ignition  of  the  end  gases.  Additionally,  developing  a  controlling 
strategy  in  HCC1  engines  (homogeneous  charge  compression  ignition) 
requires  an  understanding  of  the  auto-ignition  characteristics  of  the 
combustible  mixture  [57].  Auto-ignition  is  the  chemical  reaction  of  a 
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Fig.  8.  Ignition  temperature  as  a  function  ofXH2-  Adapted  from  Fotache  et  al.  [54]. 
Conditions:  1  atm,  strain  rate  x=150  s-1. 
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combustion  mixture,  which  releases  energy  at  a  rate  sufficient  to 
sustain  combustion  without  an  external  energy  source,  such  as  a 
spark.  Typical  auto-ignition  involves  physical  processes  (such  as  heat 
conduction,  diffusion  and  mixing  of  reactants)  and  chemical  processes 
which  have  characteristic  times  that  combine  to  form  an  ignition  delay 
time  [58],  Ideally,  ignition  delay  time  is  defined  as  the  characteristic 
chemical  time  for  ignition;  thus  experiments  are  designed  to  minimize 
the  time  associated  with  the  physical  process.  One  way  is  to  use  the 
homogeneous  mixture  in  which  diffusion  of  heat  and  mass  is  absent. 
In  this  way,  the  auto-ignition  process  is  governed  by  chemical  kinetics 
of  pre-flame  reactions.  During  this  period  of  time,  there  is  a  rapid 
depletion  of  the  primary  fuel,  high  radical  concentration  and  an 
exponential  rise  in  pressure  and  temperature;  all  these  processes  are 
governed  by  chemical  kinetics.  Thus  ignition  delay  time  is  not  only 
important  in  understanding  the  knock  phenomenon  in  ICEs  and 
developing  HCCI  controlling  strategy,  but  it  has  been  frequently  used 
to  validate  chemical  kinetic  mechanisms. 

Ignition  delay  times  of  hydrocarbons  or  hydrogen  have  been 
measured  using  a  variety  of  devices,  such  as  the  flow  reactor 
[59,60],  rapid  compression  machine  [61,62],  and  the  shock  tube 
method  [63-75].  The  shock  tube  technique  is  the  most  widely 
used  experimental  technique  for  ignition  delay  time  measurement 
because  it  can  generate  a  high  temperature  environment  under 
well-controlled  test  conditions. 

4.3.  Non-ideal  shock  conditions 

Ideal  shock  condition  assumptions  (constant  internal  energy,  U, 
and  volume,  V )  used  to  determine  the  state  of  the  gas  in  the 
reflected  shock  region  do  not  always  hold  because  of  the  growth  of 
the  boundary  layer  along  the  shock  tube.  Fig.  9(a)  shows  the 
typical  ignition  delay  time  definition  using  the  reflected  shock 


Fig.  9.  Typical  pressure  profile  (black  solid  line)  and  OH*  emission  signal  (red  dash 
line)  show  ignition  delay  in  lean  hydrogen/oxygen/argon  mixture  at  T=  1065  K  and 
p=0.95  MPa.  Constant  pressure  sustained  at  about  1700  ps  [65],  (For  interpretation 
of  the  references  to  color  in  this  figure  caption,  the  reader  is  referred  to  the  web 
version  of  this  paper.) 
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wave  method.  The  ignition  delay  time  is  defined  as  the  interval 
between  the  arrival  of  the  reflected  shock  wave  and  the  onset  of 
ignition.  The  onset  of  ignition  can  be  identified  by  the  steep 
increase  in  pressure,  as  well  as  the  OH*  emission.  It  is  seen  that 
even  for  non-reactive  argon-diluted  hydrogen/nitrogen,  there  is 
still  a  4%/ms  rise  in  the  pressure  rate.  This  pressure  increase  exists 
inherently  due  to  the  shock  tube  boundary  layer  effect  and  thus 
increases  the  temperature.  For  weak  ignition,  where  there  is  no 
steep  increase  in  pressure  (Fig.  9(b)),  the  ignition  delay  time 
should  be  defined  based  on  the  OH*  emissions  [65-68,74], 

Fig.  10  shows  the  ignition  delay  times  of  stoichiometric  hydro¬ 
gen  in  argon-diluted  oxygen  reported  in  the  literature.  Measure¬ 
ments  of  Zhang  et  al.  [65],  Pang  et  al.  [68],  and  Herzler  and 
Naumann  [76]  were  plotted  as  representative  data,  in  comparison 
with  numerical  predictions,  by  using  several  kinetic  mechanisms, 
including  GR1  3.0,  Li  et  al.’s  mechanism  [77],  USC  Mech,  and 
Konnov  et  al.’s  mechanism  [78],  It  is  seen  that  at  temperatures 
higher  than  around  1050  K,  predictions  from  different  kinetic 
mechanisms  agreed  well  with  each  other  and  also  with  experi¬ 
mental  measurements.  However,  as  the  temperature  was 
decreased,  GR1  3.0  and  Li’s  mechanism  over-predicted  the  ignition 
delay  times  while  the  USC  Mech  showed  a  lower  prediction  than 
the  measurements.  The  difference  in  predicting  the  ignition  delay 
times  among  different  kinetic  mechanisms  indicates  that  the  rate 
for  lower  temperature  reactions  requires  further  investigation. 
When  the  temperature  was  further  decreased  to  around  960  K, 
predictions  from  all  the  three  kinetic  mechanisms  were  higher 
than  the  measurements,  and  the  difference  even  exceeded  an 
order  of  magnitude.  The  deviation  between  the  predictions  of  the 
kinetic  mechanism  and  the  experimental  measurements  has  been 
explored.  Petersen  et  al.  [73]  provided  ignition  delay  time  of  syn¬ 
gas  air  mixtures  and  the  same  behavior  was  observed.  They 
attempted  to  attribute  the  discrepancy  to  the  uncertainty  in  rate 
of  certain  reactions,  such  as  C0+H02=C02+0H  in  the  tempera¬ 
ture  range.  However,  reasonable  rate  adjustments  do  not  resolve 
the  disagreements.  Dryer  and  Chaos  [66]  argued  that  the  hydro- 
gen-oxygen  chemical  induction  process  is  significantly  perturbed 
by  non-homogeneous  effects,  which  are  hard  to  remove  in 
research  experiments  in  the  regime  of  interest  (low  temperature 
and  high  pressure). 

Recently,  Pang  et  al.  [68]  measured  lower  temperature  ignition 
delay  times  of  hydrogen  in  argon-diluted  oxygen.  The  effects  of 
shock  tube  facility-dependent  gas  dynamics  and  pre-ignition 
energy  release  were  believed  to  cause  the  deviation  between 
measurements  and  kinetic  mechanism  predictions.  When  non¬ 
ideal  shock  conditions  (modified  U,  V  assumptions),  which  include 
the  measured  pressure  history,  were  used  to  determine  the  state 


Fig.  10.  Comparison  of  calculated  ignition  delays  using  current  kinetic  models  with 
ideal  U  and  V  assumption  to  that  of  the  experimental  data.  All  ignition  delays 
normalized  to  conditions  of  0.4  MPa. 


of  the  gas  behind  the  reflected  shock  wave,  the  deviations 
between  the  model  predicted  and  experimentally  measured  igni¬ 
tion  delay  times  were  finally  resolved  [68]. 

There  have  been  extensive  investigations  on  shock  tube  measure¬ 
ments  of  ignition  delay  times  of  hydrocarbons  or  hydrogen  [63-75], 
However,  significantly  less  work  has  been  conducted  for  hydrogen 
enriched  hydrocarbons.  Lifshizt  et  al.  [79]  were  among  the  first  to 
experimentally  measure  the  high  temperature  ignition  delay  times  of 
methane/hydrogen  mixtures  using  a  shock  tube.  Subsequently, 
Cheng  and  Oppenhiem  [80]  investigated  ignition  characteristics  of 
methane/hydrogen  mixtures  and  proposed  an  empirical  correlation 
with  the  ignition  delay  time  of  hydrogen  and  methane, 
t=  tch4Xh2)th”-  This  heuristic  formula  can  qualitatively  predict  the 
changes  in  activation  energy  at  high  temperature  regions.  Gersen 
et  al.  [57]  measured  the  ignition  delays  of  stoichiometric  methane/ 
hydrogen  mixtures  with  varying  hydrogen  fractions  from  0  to  100% 
at  1.5-7.0  MPa  in  a  rapid  compression  machine.  The  results  showed 
that  the  promoted  effect  of  the  hydrogen  addition  on  methane 
ignition  was  only  marginal  when  the  hydrogen  fraction  was  less 
than  20%,  while  the  ignition  delays  decreased  considerably  when  the 
hydrogen  fraction  was  over  than  50%.  Furthermore,  Gersen  et  al.  [57] 
proposed  a  mixing  relation  based  on  the  mixing  expression  (1), 
suggested  by  Cheng  and  Oppenhiem  [80],  and  it  can  approximate  the 
their  experimental  data  veiy  well.  Recently,  the  Petersen  group 
[72,81]  experimentally  measured  the  ignition  delays  of  lean 
methane/hydrogen  mixtures  with  ratios  of  80/20  and  60/40  by  using 
reflected  shock  waves.  The  significantly  reduced  ignition  delays  were 
also  observed  when  hydrogen  was  added.  Furthermore,  they  found 
that  the  addition  of  hydrogen  did  not  seem  to  shift  the  dominant 
kinetic  in  their  study. 

Huang  et  al.  [82]  conducted  both  experimental  and  numerical 
studies  on  auto-ignition  and  chemical  kinetics  of  stoichiometric 
methane/hydrogen/air  mixtures  with  varying  hydrogen  fractions 
from  15%  to  35%  at  high  pressure  (1.6  MPa  and  4.0  MPa)  in  a  shock 
tube.  They  reported  that  the  effect  of  the  hydrogen  addition  on 
ignition  methane  was  presented  and  the  reduced  extent  of  ignition 
delays  decreased  with  the  decrease  in  temperature.  Chaumeix 
et  al.  [83]  investigated  the  shock  tube  ignition  delay  and  detona¬ 
tion  characteristic  of  methane/hydrogen  mixtures;  they  concluded 
that  the  presence  of  methane  significantly  inhibits  the  detonation 
of  the  combustible  mixtures  due  to  lowered  ignition  temperature 
and  longer  ignition  delay  of  methane. 

Herzler  and  Naumann  [76]  conducted  the  experimental  and 
numerical  study  on  ignition  delays  of  methane/ethane/hydrogen 
mixtures  with  varying  hydrogen  fractions  from  0  to  100%  in  a 
shock  tube.  Their  results  showed  that  the  current  mechanisms 
could  not  well  represent  the  reduction  of  global  activation  energy 
at  relative  low  temperature  for  high  fractions  of  hydrogen  in  the 
fuel  mixtures. 

More  recently,  a  more  systematic  investigation  on  auto-ignition 
characteristics  of  methane/hydrogen  mixtures  was  conducted  by 
Zhang  et  al.  [65],  and  three  auto-ignition  regimes  were  identified, 
depending  on  the  hydrogen  fraction:  the  methane  chemistry  domi¬ 
nated  ignition  (MCD1),  when  the  hydrogen  fraction  was  less  than 
40%;  the  combined  chemistry  of  methane  and  hydrogen  dominated 
ignition  (CCMHDI),  when  the  mole  fraction  of  hydrogen  was  about 
60%;  and  the  hydrogen  chemistiy  dominated  ignition  (HCDI),  when 
the  hydrogen  fraction  was  more  than  80%.  Detailed  discussion  of 
auto-ignition  characteristics,  and  the  corresponding  reaction  kinetics 
in  each  regime,  based  on  the  work  by  Zhang  et  al.  follows. 

4.2.  Reaction  system  in  MCDI  regime 

Zhang  et  al.  [65]  found  that  for  hydrogen  fractions  XH2  with  a  less 
than  40%  methane/hydrogen  mixture,  the  ignition  characteristics  of 
methane/hydrogen  resembles  that  of  pure  methane  and  the  ignition 
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process  is  governed  by  the  kinetics  of  methane;  thus  this  regime  is 
defined  as  the  methane-chemistry-dominated-ignition  (MCDI) 
regime.  In  this  regime,  the  ignition  delay  time  dependence  on 
pressure  and  equivalence  ratio  is  specified  in  the  following. 

4.2.1.  Effect  of  pressure  in  MCDI  regime 

Fig.  11  shows  the  ignition  delay  time  of  methane/hydrogen  in 
Ar  diluted  oxygen  forXH2  =  0,  20%,  and  40%,  an  equivalence  ratio  of 
0.5,  and  at  elevated  pressures.  For  pure  methane  fuel,  as  shown  in 
Fig.  11(a),  the  logarithmic  ignition  delay  time  is  a  linear  function  of 
the  inverse  temperature;  thus  an  Arrhenius-type  correlation, 
through  a  multi-regression  method,  was  formulated,  as  shown  in 
Eq.  (6),  where  R=  1.986  x  10~3  kcal/(mole  K)  is  the  universal  gas 
constant,  A,  B,  C,  and  D  are  correlation  parameters  and  Ea  is  the 
global  activation  energy  in  kcal/mole.  As  the  pressure  increases 
from  0.5  MPa  to  2  MPa,  the  global  activation  energy  is  decreased 
from  48.4  to  42.8  kcal/mole 

Tign  =Ap-  VC*02D  exp (Ea/KT)  (6) 


10000  K/  T 


10000  K/  T 

Fig.  11.  Effect  of  pressure  on  ignition  for  XH2  <  40%  in  lean  methane/hydrogen 
mixtures  with  dilution  Ar.  Adapted  from  Ref.  [65], 


Generally,  the  pressure  exponent,  B,  in  Eq.  (6),  is  a  negative 
value  of  around  0.5  for  typical  hydrocarbon  fuels,  which  implies 
that  ignition  delay  time  is  reduced  with  increasing  pressure. 
Additionally,  typical  values  of  the  global  activation  energy  of  pure 
methane  at  the  temperature  range  between  1315  and  2000  K  are 
within  the  range  of  46-54  kcal/mol  [72,75,79,80,84,85],  as  shown 
in  Table  2.  Relatively  small  activation  energy  (18  kcal/mol)  was 
reported  in  the  study  of  Huang  et  al.  [86],  It  is  noted  that  the 
dominant  chain-branching  reactions  are  quite  different  in  different 
temperature  ranges.  The  chain-branching  reactions  CH3+02*->0+ 
CH30  and  H02  +  CH3  <->  OH  +  CH30  are  more  important  when  the 
temperature  is  higher  than  1400  K.  However,  the  reactions 
CH3  +  CH302^CH30  +  CH30  and  H202  (  +  M)^OH  +  OH  (  +  M) 
dominate  the  ignition  chemistry  when  the  temperature  is  lower 
than  1100  I<  [69],  Thus  the  significant  difference  in  global  activa¬ 
tion  energy  reported  by  Huang  et  al.  is  caused  by  the  temperature 
range  (1000-1350  K),  which  they  studied.  Furthermore,  Huang 
et  al.  [86]  found  a  reversed  “S”  shaped  characteristic  in  the 
methane/air  mixture.  Petersen  et  al.  [71]  suggested  that  the 
elevated  temperature  would  not  shift  the  kinetic  regimes  for  the 
lean  methane-oxygen  mixtures.  Note  that  the  effect  of  elevated 
pressure  on  the  ignition  kinetic  was  weakened  with  an  increase  in 
the  hydrogen  mole  fraction  in  the  fuel  mixture  for  MCDI  at  the 
relative  low  temperature  region.  When  the  hydrogen  was  added 
into  methane,  as  shown  in  Fig.  11(b)  and  (c),  there  was  a  slight 
change  in  the  global  activation  energy  for  the  pressure  of  1  and 
2  MPa.  This  is  because  the  hydrogen  addition  weakened  the 
dominance  of  the  methane  chemistry.  The  similar  transition  in 
global  activation  energy  was  also  observed  by  Petersen  et  al.  [72], 
for  XH2=20%  CH4/H2  fuel  blend  at  a  pressure  of  2.1  MPa.  A  high 
temperature  activation  energy  (EaH)  of  41.4  kcal/mol  and  a  low 
temperature  activation  energy  (£aL),  of  31.1  kcal/mol  were  reported 
in  their  study.  The  transition  was  attributed  to  the  effect  of  the 
hydrogen  addition  on  the  chain-branching  during  the  ignition 
process. 

Herzler  and  Naumann  [76]  measured  the  ignition  delay 
times  of  stoichiometric  hydrogen/reference  gas  (mixtures  of  92% 


Table  2 

Global  activation  energy  of  typical  hydrocarbon  fuels. 


Source 

Experimental  conditions 

Activation 
energy  (kcal/mol) 

Zhang  et  al.  [65] 

I5=  1290-2000  K 
p5=0.5-2.0  MPa 
<f> =0.5 

Ar=80% 

42.8-48.4 

Grillo  and  Slack  [84] 

T5=  1400-1850  K 
p5=0.4  MPa 
(f>= 2.0 

Ar= 79-96% 

53.4 

Seery  and  Bowman  [85] 

T5=  1350-1900  I< 
p5= 0.15-0.4  MPa 

0= 0.2-5.0 

Ar= 53.4-78.4% 

52.3 

Petersen  et  al.  [72,75] 

T5=  1410-2040  K 
p5= 0.9-48  MPa 

0=O.5-4.O 

Ar= 89-99%,  N2  =  97.66% 

51.8 

Lifshitz  et  al.  [79] 

T5=  1500-2100  K 
p5=0.2-1.0  MPa 
(f)= 0.5-2.0 

Ar= 89-97% 

51.4 

Cheng  and  Oppenheim  [80] 

T5= 800-2400  K 
p5=0.1-0.3  MPa 

0=0.5-15 

Ar=90% 

46.4 
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methane  and  8%  ethane)  in  Argon-diluted  oxygen  (dilution  ratio 
1:5),  and  their  measurements,  along  with  numerical  predictions, 
using  USC  Mech  II  at  the  exact  conditions  presented  in  Fig.  12. 
Excellent  agreement  between  the  measurements  and  the  numer¬ 
ical  predictions  was  observed  in  the  temperature  range  studied  in 
Ref.  [76].  The  ignition  delay  time  of  the  reference  gas  (shown  in 
Fig.  12(a))  decreased  with  increasing  pressure  with  a  factor  of  p0  5, 
which  is  typical  for  hydrocarbons.  Their  data  for  the  equivalence 
ratio  of  0.5  were  not  plotted  here,  and  all  the  data  for  the  reference 
gas  can  be  fitted  by 

zigJps  =  10-  2  75  ±  0 13(p/bar)  ~  °-51  ±  0  02  /  59  ±  006 

xexp(20, 450  ±  442  K/T)  (7) 

yielding  a  global  activation  energy  of  40.46  kcal/mole. 

4.2.2.  Effects  of  equivalence  ratio  in  MCD1  regime 

Fig.  13  shows  the  ignition  delay  time  of  methane/hydrogen  in 
Ar  diluted  oxygen  for  the  equivalence  ratio  of  0.5, 1.0  and  2.0  at  the 
pressure  of  2  MPa.  It  is  seen  that  two  temperature  regions  were 
classified,  based  on  the  changes  in  global  activation  energy.  At  the 
high  temperature  region,  ignition  delay  times  increased  with  the 
increasing  equivalence  ratio,  indicating  that  the  lean  mixture  had 
higher  reactivity  during  the  ignition  process.  The  ignition  kinetic 
depends  more  on  oxygen  concentrations  at  high  temperature 
regions  because  the  chain-branching  reaction  H+02=0H-|-0  is 
facilitated  for  leaner  mixtures  (higher  oxygen  concentration),  in 
the  high  temperature  range,  and  thus  promotes  the  ignition  of  the 
mixtures  and  reduces  the  ignition  delay  times. 

When  the  temperature  is  decreased  (demonstrated  by  the  open 
symbols  in  Fig.  13),  a  similar  dependence  on  equivalence  ratio  is 
presented:  the  ignition  delays  of  mixtures  increased  with  the 
increasing  equivalence  ratio.  However,  the  overall  slope  of  the 


Fig.  12.  Effect  of  pressure  on  ignition  for  XH2  <  40%  in  stoichiometric  methane/ 
ethane/hydrogen  mixtures  with  dilution  Ar.  Experiments:  symbols:  obtained  from 
Herzler  et  al.  [76].  Simulations:  lines:  USC  2.0  Mech  [134]  with  dp/dt=4%/ms. 


Fig.  13.  Effect  of  equivalence  ratio  on  ignition  for  MCDI  mixtures  with  dilution  Ar. 
Adapted  from  Ref.  [65], 

logarithmic  ignition  delay  time  -  or  the  global  activation  energy  - 
was  significantly  lower  than  at  high  temperatures.  When  the 
temperature  was  decreased,  the  fuel-specific  reactions  became 
more  important  in  governing  the  ignition.  Specifically,  the  con¬ 
centration  of  methane  increased  with  the  increased  equivalence 
ratio,  leading  to  the  increased  production  of  methyl  radicals,  which 
in  turn  promoted  the  chain  termination  reaction  CH3  +  CH3  (  +  M) 
<->C2H6  (  +  M),  and  increased  the  ignition  delay  time.  Additionally, 
the  chain-branching  reaction  H  +  02<->0H+0  at  high  tempera¬ 
tures  has  a  higher  activation  energy  than  the  methyl  radical 
production  reactions  by  H  abstraction  from  methane,  resulting  in 
the  observed,  more  gentle  slope  at  lower  temperatures. 

4.3.  Reaction  system  in  CCMHDI  regime 

4.3.1.  Effects  of  pressure  in  CCMHDI  regime 

For  XH2  around  60%,  Zhang  et  al.  proposed  that  the  ignition 
delay  time  of  the  methane/hydrogen  in  Ar  diluted  oxygen  is 
governed  by  the  combined  chemistry  of  methane  and  hydrogen 
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Fig.  14.  Effect  of  pressure  on  ignition  delay  for  CCMHDI  mixture.  Adapted  from  Ref.  [65]. 


Fig.  15.  Effect  of  equivalence  ratio  on  ignition  delay  for  CCMHDI  mixture.  Adapted 
from  Ref.  [65], 

dominated  ignition  (CCMHDI).  Fig.  14  shows  the  ignition  delay 
time  of  methane/hydrogen  mixtures  for  the  equivalence  ratio  of 
0.5  and  a  pressure  of  0.5,  1  and  2  MPa.  Unlike  the  pressure 
behavior  in  the  MCDI  shown  in  Fig.  11,  the  ignition  delay  times 
at  the  three  pressures  investigated  approach  the  same  value, 
which  implies  that  the  pressure  imposes  only  subtle  influence 
on  the  ignition  delay  time.  The  pressure  dependence  is  different 
from  those  of  MCDI  mixtures.  An  Arrhenius-type  expression  with 
R2  =  0.97  is  correlated  using  the  multiple  linear  regression  method 
based  on  the  normalized  pressure  of  1.0  MPa 

rign//rs=  10-41  ±023(p/bar)-000629  exp(19,100  +  623  K/T)  (8) 

The  insensitivity  of  the  ignition  delay  time  on  pressure  for 
CCMHDI  mixtures  was  reasonable  because  for  methane  the 
pressure  imposed  a  promoting  effect  on  ignition,  as  seen  in 
Eq.  (7).  However,  when  the  pressure  was  increased  for  argon- 
diluted  hydrogen/oxygen  mixtures,  the  termination  reaction 
H+02=H02  is  favored.  Thus  the  net  effect  was  the  observed  weak 
dependence  of  the  ignition  delay  time  on  pressure. 

4.3.2.  Effects  of  equivalence  ratio  in  CCMHDI  regime 

For  the  CCMHDI  mixture,  Fig.  15  shows  ignition  delay  time  for 
the  equivalence  ratio  of  0.5,  1.0  and  2.0.  In  the  high  temperature 
region,  the  ignition  delay  time  for  tp  =1.0  is  slightly  higher  than 
<p= 0.5,  while  the  ignition  delay  for  (p= 2.0  is  significantly  higher. 
This  non-linear  variation  of  the  logarithmic  ignition  delay  time 
with  the  equivalence  ratio  indicates  that  similar  correlation  as 
Eq.  (7)  for  the  ignition  delay  time  as  a  function  of  the  equivalence 
ratio  is  not  available. 

In  the  middle-low  temperature  region,  ignition  delay  times 
are  insensitive  to  the  change  of  equivalence  ratio.  A  similar 


phenomenon  was  reported  by  Gersen  et  al.  [57];  the  reason  is 
that  ignition  is  now  governed  by  the  combined  chemistry  of 
methane  and  hydrogen.  The  typical  reactions  that  favor  ignition 
are  the  chain-branching  reaction  H+02<-»0H+0,  while  the  most 
effective  ignition  inhibiting  reaction  is  H  radical  consumption 
reaction  CH4+H<->CH3  +  H2.  As  the  equivalence  ratio  is  increased, 
the  hydrogen  concentration  increases,  leading  to  increased  pro¬ 
duction  of  H  radicals  and  consequently,  decreased  ignition  delay 
time.  On  the  other  hand,  the  concentration  of  methane  also 
increases,  which  consumes  the  H  radical  and  increases  the  ignition 
delay  time.  Thus  the  insensitivity  of  delayed  ignition  time  on  the 
equivalence  ratio  in  the  temperature  range  is  caused  by  the 
combined  chemistry  of  methane  and  hydrogen. 

4.4.  Reaction  system  in  HCDI  regime 

ForXH2  higher  than  80%,  it  is  believed  that  the  ignition  behavior 
of  methane/hydrogen  fuel  mixtures  is  dominated  by  hydrogen 
chemistry,  and  Zhang  et  al.  [65]  proposed  that  ignition  behavior  of 
this  type  of  fuel  mixtures  to  be  hydrogen  chemistry  dominated 
ignition  (HCDI). 

4.4.1.  Effect  of  pressure  in  HCDI  regime 

Fig.  16  shows  the  ignition  delay  time  for  XH2=80%  methane/ 
hydrogen  fuel  mixtures  in  argon-diluted  oxygen  at  the  equivalence 
ratio  of  0.5,  and  elevated  pressures.  It  is  seen  that  when  ignition 
delay  time  versus  inverse  temperature,  there  is  a  significant 
transition  in  slope,  indicating  a  change  in  the  overall  activation 
energy.  Skinner  and  Ringrose  [87]  investigated  the  ignition  delay 
time  of  pure  hydrogen  in  oxygen  (argon-diluted);  their  results  also 
showed  two  transitions  in  overall  activation  energy  with  increased 
temperature.  At  higher  temperatures  (T >  1200  K),  the  ignition 
delay  time  exhibited  a  negative  dependence  on  pressure,  indicating 
that  the  overall  reactivity  increased  with  the  increasing  pressure.  An 
opposite  pressure  dependence  is  presented  at  intermediate  tempera¬ 
ture.  This  observation  is  consistent  with  the  study  of  Herzler  and 
Naumann  [76]  for  methane/ethane/hydrogen  mixtures.  The  reason  is 
that  when  hydrogen  fraction  is  over  80%,  ignition  is  governed 
primarily  by  hydrogen  chemistry.  The  competition  between  reactions 
H+02  (  +  M)<->H02  (  +  M)  and  H+02«OH  +0  is  believed  to  result 
in  complex  pressure  dependence.  Specifically,  the  ignition  delay  time 
of  the  HCDI  system  at  2.0  MPa  and  for  T <  1200  I<  is  higher  than  at 
0.5  and  1.0  MPa  because  reaction  H+02  (  +  M)<->H02  (  +  M)  is 
facilitated  at  higher  pressures.  While  for  higher  temperatures  and 
typical  hydrocarbon  chemistry  dominated  ignition,  the  higher  abso¬ 
lute  concentrations  due  to  higher  pressures  will  promote  the  chain¬ 
branching  reaction  H+02<->0H+0,  thus  decreases  ignition  delay 
time,  as  observed  in  Fig.  16  for  temperatures  higher  than  1100  K  and 
in  Figs.  11-14  for  typical  hydrocarbon  chemistry  dominated  ignition. 


Fig.  16.  Effect  of  pressure  on  ignition  delay  for  HCDI  mixtures  (20%  methane/80% 
hydrogen).  Adapted  from  Ref.  [65], 
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In  addition,  Fig.  16  shows  that  the  crossing  of  data  at  1.0  and  0.5  MPa 
shifted  to  even  lower  temperatures  (1100  K).  This  observation  can,  by 
the  same  token,  be  explained  by  the  fact  that  reaction  H+02  (  +  M) 
<-*H02  (  +  M)  -  if  preferably  facilitated  over  the  reaction  H+02<-> 
OH+O  at  relatively  higher  pressure  and  lower  temperature. 

4.4.2.  Effect  of  equivalence  ratio  in  HCDI  regime 

Fig.  17  shows  the  ignition  delay  time  for  80%  reference  gas/20% 
hydrogen  fuel  mixtures  in  argon-diluted  oxygen  at  0.1  MPa  and 
the  equivalence  ratio  of  0.5  and  1.0  [76],  The  ignition  delay  time 
increases  with  the  increase  of  equivalence  ratio  because,  in  this 
temperature  range,  reducing  the  oxygen  concentration  retards  the 
chain-branching  reaction  H  +  02<-»0H  +  0  and  consequently  inhi¬ 
bits  ignition. 

However,  Herzler  and  Naumann  [76]  showed  that  at  higher 
pressure  (1.6  MPa)  a  weak  influence  of  equivalence  ratio  occurred 


Fig.  17.  Effect  of  equivalence  ratio  on  ignition  delay  for  HCDI  mixture  at  low 
pressure.  Experiments:  symbols:  obtained  from  Herzler  et  al.  [76].  Simulations: 
lines:  USC  2.0  Mech  [134]  with  dp/dt=4%/ms. 


Fig.  18.  Effect  of  equivalence  ratio  on  ignition  delay  for  HCDI  mixtures.  Adapted 
from  Ref.  [65], 


in  the  ignition  delay  time  for  the  HCDI  mixture,  as  shown  in 
Fig.  18(a).  In  the  lower  temperature  region  (T <  1160  K),  ignition 
delay  times  at  the  equivalence  ratio  of  0.5  and  1.0  were  almost 
identical.  When  the  temperature  was  increased  to  1250  K,  the 
ignition  delay  time  at  <p= 0.5  was  slightly  lower  than  at  <p=  1.0. 
This  observation  is  consistent  with  the  results  of  Zhang  et  al.  [65], 
who  measured  ignition  delay  times  for  XH2  =  80%  hydrogen/ 
methane  in  argon-diluted  oxygen  at  1.8  MPa  for  cp= 0.5,  1.0,  and 
2.0,  as  shown  in  Fig.  18(b).  It  is  seen  that  in  lower  temperature 
regions,  the  ignition  delay  time  for  three  equivalence  ratios  was 
almost  the  same,  while  at  higher  temperatures,  the  ignition  delay 
time  increased  with  the  increased  equivalence  ratio. 


5.  Laminar  flame  speed  and  flame  front  instability 

5.3.  Laminar  flame  speed 

Laminar  flame  speed  is  a  physiochemical  property  of  a  com¬ 
bustible  mixture,  which  characterizes  its  diffusivity,  reactivity  and 
exothermicity  [25].  Accurately  measured  laminar  flame  speed  can 
be  used  to  validate  chemical  kinetic  mechanisms  [88-90]  in 
engine  design  [91]  and  in  turbulent  combustion  modeling  [92-94], 

Hydrocarbon  oxidation  mechanisms  have  been  reported  and 
reviewed  extensively  [89,90];  and  they  have  been  greatly 
improved,  enabling  the  numerical  prediction  of  laminar  flame 
speed  of  each  single-fuel/oxidizer  system.  However,  there  is  still 
no  method  for  calculating  the  laminar  flame  speed  of  hydrogen/ 
hydrocarbon/air  mixtures,  because  the  oxidation  mechanisms  of 
fuel  blends  cannot  be  the  simple  combinations  of  the  mechanisms 
of  all  the  constituent  fuels.  Thus  experiments  are  required  to 
obtain  the  laminar  flame  speed  data  for  hydrogen/hydrocarbon 
fuel  blends. 

Methods  for  measuring  laminar  flame  speed  include  the 
stagnation  plane  flame  method  [95],  the  heat  flux  method 
[96,97],  and  the  combustion  bomb  method  [98-103],  The  stagna¬ 
tion  plane  flame  method  utilizes  the  stable  flow  of  premixed 
combustible  mixtures  and  the  flow  velocity  was  profiled  through 
which  the  strain  rate  is  extracted.  When  the  flow  velocity  was 
extrapolated  to  zero  strain  rate,  the  laminar  flame  speed  was 
determined.  The  heat  flux  method  extrapolates  the  heat  flux  to 
zero  to  obtain  the  adiabatic  flame  speed.  While  these  two  methods 
are  easy,  because  the  parameters  to  be  measured  are  stable 
quantities.  However,  typically  the  flame  speed  at  high  pressures 
is  difficult  to  be  measured  through  the  stagnation  plane  flame 
method  or  the  heat  flux  method  because  at  high  pressures,  stable 
combustion  is  inherently  difficult  to  achieve.  The  combustion 
bomb  method  introduces  outwardly  propagating  spherical  flame 
and  has  been  widely  used  because  of  its  simple  flame  configura¬ 
tion,  well-defined  flame  stretch  rate  and  well-controlled  experi¬ 
mentation  [98,99],  making  it  the  most  widely  used  experimental 
technique  for  flame  speed  measurements.  The  first  outwardly 
propagating  spherical  flame  was  imaged  and  reported  by  Manton 
et  al.  [104];  however,  their  focus  was  on  flame  morphology,  e.g., 
flame  cellularity  caused  by  preferential  diffusion.  The  first  flame 
speed  measurements  with  schlieren  observation  were  reported  in 
1969  by  Palm-Leis  and  Roger  105];  they  attributed  the  changes  in 
flame  speed  in  spherical  explosions  to  flame  stretch.  In  addition, 
most  previous  experimental  reports  on  laminar  flame  speeds 
measurements  used  hydrogen/air  [99,106-110],  or  hydrocarbons/ 
air  [96-98,111-125],  mixtures  and  there  is  relatively  less  literature 
on  the  effects  of  hydrogen  addition  on  the  laminar  flame  speed  of 
hydrocarbon/air  mixtures  [95,100-103,126-128].  Although  the 
definitions  of  hydrogen  addition  extent  used  in  these  studies  were 
different,  qualitatively,  all  these  studies  show  that  laminar  flame 
speed  increases  when  more  hydrogen  is  added.  The  definition  of 
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each  hydrogen  addition  parameter  and  presentations  of  the  flame 
speed  variation  of  each  parameter  follows. 


5.1.3.  Laminar  flame  speed  versus  XH2 

Currently,  most  researchers  use  XH2  to  quantify  the  hydrogen 
addition  extent  because  is  the  simplest  quantification  parameter 
[100-102,127,129-132],  possessing  significant  merit  in  engineer¬ 
ing  applications,  such  as  in  hydrogen/natural  gas  spark  ignition 
engines  [133], 


5.1.3. 3.  H2/natural  gas  or  methane.  Huang  et  al.  [100]  studied  the 
effects  of  hydrogen  addition  on  laminar  flame  speeds  of  natural 
gas  at  1  atm  and  room  temperature;  they  found  that  normalized 
flame  speed  increased  quasi-exponentially  with  XH2-  An  empirical 
correlation  between  the  laminar  flame  speed  and  XH2  was 
deduced  from  the  experimental  observations,  as  shown  in  the 
following  equation: 


S°0/->,Xh2)-5°(</>,0) 
S°(4>,1)  — S°(4>,0) 


000737  “Kora)*0-00334 


O) 


Methane  is  a  major  component  of  natural  gas  and  the  laminar 
flame  speed  of  hydrogen/methane/air  mixtures  has  been  compu¬ 
tationally  [126],  and  experimentally  [127]  studied.  Three  distinct 
regimes  of  flame  speed  behavior,  with  increasing  XH2  was  identi¬ 
fied,  as  shown  in  Fig.  19.  When  0<XH2<0.6  (Regime  1),  the 
laminar  flame  speed  increased  slightly  with  increasing  XH2,  indi¬ 
cating  that  in  this  regime,  the  combustion  was  dominated  by 
methane.  When  XH2  is  higher  than  0.8  (Regime  III),  the  laminar 
flame  speed  decreased  dramatically  with  decreasing  XH2,  indicat¬ 
ing  that  the  presence  of  small  amount  of  methane  can  greatly 
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Fig.  19.  Laminar  flame  speed  of  hydrogen/methane/air  mixtures  as  a  function  of 
XH2  at  ^=1.0.  Scatters:  measurements  (data  source:  Ref.  [127]),  lines:  calculations 
(data  source:  Ref.  [126]). 


inhibit  hydrogen/air  flame  propagation.  When  0.6<XH2<0.9 
(Regime  II),  the  flame  speed  increased  by  increasing  XH2  is 
respectively  more  and  less  significant  than  Regimes  I  and  III. 

5.3. 3.2.  H^propane.  The  laminar  flame  speed  of  stoichiometric 
hydrogen/propane/air  mixtures  was  experimentally  reported  by 
Milton  and  Keck  [91  ],  using  the  closed  bomb  method.  Their  results 
were  obtained  through  thermodynamic  calculation  based  on  the 
pressure  history;  their  calculated  flame  speed  has  complex 
dependence  on  the  initial  temperature  and  pressure,  as  shown 
in  Fig.  20.  When  the  H/C  ratio  is  8  (equivalent  to  XH2=0.89), 
laminar  flame  speed  is  much  higher  than  when  H/C  ratio  equals  4 
(equivalent  to  XH2  =  0.67). 

Recently,  Law  and  Kwon  [128]  studied  the  hydrocarbon  sub¬ 
stitution  effect  on  atmospheric  hydrogen  air  propagation  with  an 
outwardly  spherical  flame.  They  measured  the  flame  speed  of 
atmospheric  hydrogen/air  with  small  amounts  of  hydrocarbon 
substitution;  the  laminar  flame  speeds  were  plotted  against  the 
specifically  defined  SFHC  (stoichiometric  fraction  of  hydrocarbon, 
defined  as  5XaH8/(0.5XH2+5Xc3H8).  Ref-  [128]).  Tang  et  al.  also 
measured  the  flame  speed  of  propane  with  an  increasing  hydrogen 
addition  extent  at  atmospheric  pressure,  room  temperature  [102], 
and  elevated  pressures  and  temperatures  [101],  The  objectives  of 
Law  and  Kwon  and  Tang  et  al.  were  quite  different:  Law  and  Kwon 
[128]  focused  on  the  hydrocarbon  addition  effect  on  hydrogen/air 
flame  propagation,  while  Tang  et  al.  [101]  -  with  a  practical  and 
simply-defined  hydrogen  addition  extent  -  measured  the  flame 
speed  of  propane  with  a  hydrogen  addition.  The  SFHC  defined  in 
Ref.  [128]  were  converted  to  XH2  in  Ref.  [101],  and  their  results  of 
propane/hydrogen  flame  speeds  were  plotted  in  Fig.  21  as  a 
function  of  XH2.  It  clearly  shows  that  for  XH2  larger  than  80%, 
laminar  flame  speed  increased  dramatically  with  the  increase  of 
XH2-  The  results  of  Law  and  Kwon  covered  the  high  XH2  (low  SFHC), 
range  and  the  results  of  Tang  et  al.  [101],  being  quite  the  opposite, 
covered  the  whole  range,  but  very  sparsely  at  high  XH2-  However,  it 
was  found  that  a  simple  fitting  of  the  results  of  Tang  et  al.  overlaps 
the  results  of  Law  and  Kwon,  thus  unifying  the  definitions  of  XH2 
and  SFHC. 

The  results  of  Tang  et  al.  [101]  were  then  compared  with  the 
computationally  predicted  values.  The  mechanism  they  used  was  the 
recently  developed  and  highly  cited  USC  Mech  II  [134],  for  high 
temperature  oxidation  of  H2/CO/Cl-C4  compounds.  It  showed  that 
agreement  between  the  experimental  and  computational  results  was 
excellent  at  low  atmospheric  pressures,  while  discrepancies  still 
existed  for  higher  pressures.  This  shows  the  irreconcilable  difference 
between  experimental  measurements  and  computational  predictions. 
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Fig.  20.  Laminar  flame  speed  of  stoichiometric  hydrogen/propane/air  mixtures  as  a 
function  of  temperature  and  pressure.  Data  source:  Ref.  [91  ]. 
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Fig.  21.  Flame  speed  of  Law  and  Kwon  [128]  and  Tang  et  al.  [101]  as  a  function  of  XH2- 


Fig.  22.  Experimental  (points)  and  calculated  (line)  laminar  flame  speeds  as  a 
function  of  5  atm  and  overall  equivalence  ratio  of  1.0,  adapted  from  Law  et  al.  [135J. 


The  results  of  Law  et  al.  [135]  were  also  presented  here,  in  Fig.  22.  The 
figure  clearly  shows  that  both  experimentally  measured  and  compu¬ 
tationally  predicted  flame  speed  decreased  with  the  increase  of 
propane  addition.  At  a  lower  mixing  ratio  (higher  hydrogen  fraction), 
the  agreement  between  the  experiments  and  the  computations 
becomes  bad  because  the  mechanism  of  Qin  et  al.  [136]  was  devel¬ 
oped  for  C1-C2  hydrocarbons  and  could  not  degenerate  to  the 
hydrogen  oxidization  mechanism  as  the  mixing  ratio  decreases 
(hydrogen  fraction  increases). 

The  dependence  of  laminar  flame  speed  on  temperature  and 
pressure  is  discussed  in  Ref.  [101].  Laminar  flame  speed  decreases 
monotonically  with  the  increase  of  Pu  and  increases  quasi-linearly 
with  the  increase  of  Tu.  As  the  pressure  increases,  the  intensity 
of  the  temperature-sensitive,  two-body,  branching  reaction:  H  + 
02  *-♦011  —  0  is  approximately  fixed,  due  to  the  insensitivity  of 
adiabatic  temperature  to  the  increase  in  pressure,  while  the  three- 
body,  temperature-insensitive,  inhibiting  reaction:  H+02+M<-> 
H02  +  M  is  enhanced,  and  a  retarding  effect  is  therefore  imposed 
on  the  overall  progress  of  the  reaction  with  increasing  pressure. 
However,  this  effect  is  not  great  enough  to  be  responsible  for 
the  trend  of  the  flame  speed  to  decrease  with  the  increasing 
pressure  because  the  pressure  effect  of  the  chain  mechanism  can 
be  quantified  through  the  overall  reaction  order  n  (typically  in 
the  range  of  1-2),  which  can  be  locally  defined  as  /°~pn/2.  Thus 
the  flame  speed  ~/°//?u~pn/2~1  shows  that  the  decrease  in  flame 
speed  with  the  increasing  pressure  is  actually  caused  by  increasing 
pressure  [137-139].  The  increase  in  the  upstream  temperature 
increases  the  adiabatic  temperature,  which  influences  the  reaction 


rate;  the  dependence  is  more  sensitive  for  larger  values  of  Tu 
because  of  the  Arrhenius  factor. 


5.3. 1.3.  H2/higher  hydrocarbons.  Mandilas  et  al.  [140]  measured  the 
laminar  flame  speed  of  hydrogen/iso-octane/air  mixtures.  About 
5%  by  mass  of  hydrogen  was  added  (equivalent  to  XH2  =  0.75)  into 
the  base  fuel  and  it  was  observed  that  the  laminar  flame  speed  of 
hydrogen/iso-octane/air  mixtures  was  higher  than  those  of  pure 
iso-octane/air  mixtures  over  the  equivalence  ratio  range  between 
0.8  and  2.0. 


5.3.2.  Laminar  flame  speed  versus  RH 

5.12.1.  Definition  of  RH.  In  1984,  Yu  et  al.  [95]  used  the  symmetric 
counter-flow  flame  to  accurately  measure  the  flame  speeds  of 
propane-air  with  stoichiometrically  small  amounts  of  hydrogen 
addition.  In  their  study,  it  was  assumed  that  the  hydrogen  only 
presents  in  small  amounts  and  is  completely  consumed  by  part  of 
the  total  air;  in  consequence  the  remaining  air  is  the  oxidizer 
defining  the  mixture  stoichiometry.  Thus,  the  stoichiometry  was 
defined  through  the  mole  fraction  of  hydrogen  (CH),  propane  (CF), 
and  air  (CA)  as 


CF/[CA-CH/(CH/CA)stoic] 

(Cp/CA  )st0ic 


(10) 


where  the  subscript  represents  the  stoichiometric  value  (for  instance, 
(CH/CA)stoic  is  the  stoichiometric  fuel-to-air  molar  ratio);  the  remain¬ 
ing  air  for  hydrocarbon  oxidation  is  represented  by  the  term  in 
brackets. 

The  extent  of  hydrogen  addition  was  then  correspondingly 
defined  as 


CH  +  CH/(CH/CA)stoic 
Cf  +  [CA  —  Ch/ (CH/CA)stoic] 


(11) 


It  should  be  pointed  out  that  the  parameters  (f> F  and  RH  indicate 
the  mixture  concentration  and  the  relative  amount  of  added 
hydrogen.  In  addition,  it  is  easy  to  discover  that  the  two  para¬ 
meters  do  not  represent  the  actual  stoichiometry  during  the 
reaction,  but  they  do  facilitate  data  reduction  and  correlation, 
which  will  be  discussed  in  detail  in  the  following. 

5.12.2.  Linearity  between  the  laminar  flame  speed  and  RH.  Fig.  23 
shows  the  laminar  flame  speed  of  hydrogen  enriched  methane  and 
propane-air  mixtures  as  a  function  of  RH  at  1  atm  and  room 
temperature  [95],  The  scatters  were  the  experimental  results  and 
it  is  seen  that  the  linear  fitting  (lines)  correlates  the  laminar  flame 
speed  very  well  with  RH.  Yu  et  al.  [95]  derived  the  following 
correlation  for  the  laminar  flame  speed  of  hydrogen  enriched 
methane  and  propane  fuel  mixtures: 

S0u((pF,RH)  =  S°u((PF,0)+83RH  (12) 


Eq.  (12)  indicates  that  the  laminar  flame  speed  is  a  linear 
function  of  RH;  in  addition,  the  slope  of  linearity  is  independent  of 
the  effective  equivalence  ratio,  which  means  that  the  hydrogen 
addition  effect  is  quantitatively  the  same  for  different  effective 
equivalence  ratios.  This  correlation  has  elicited  substantial  follow¬ 
up  investigations  because  of  the  interest  of  hydrogen  combustion 
and  its  simplicity. 

Sher  and  Ozdor  [141]  have  determined  the  laminar  flame  speed 
of  hydrogen  enriched  butane/air  mixtures  by  thermodynamic 
calculations  based  on  the  flow  rate  of  the  mixture  and  heat  carried 
away  by  the  cooling  water. 

Since  the  mid-1980s,  the  accuracy  of  the  determination  of 
laminar  flame  speed  has  enabled  important  improvements  due  to 
the  development  of  experimental  techniques  such  as  high  speed 
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Fig.  23.  Laminar  flame  speed  of  hydrogen  enriched  (a)  methane  and  (b)  propane  as 
a  function  of  RH.  Data  re-plotted  from  Yu  et  al.  [95]. 

photography,  more  sophisticated  data  processing  methods  such  as 
non-linear  extrapolation  techniques  for  removal  of  stretch  effect, 
and  the  routinely  used  computational  simulations  of  ID  freely 
propagating  planar  flame.  Using  the  above  concepts,  Tang  et  al. 
[103]  experimentally  re-examined  and  extended  the  linear  corre¬ 
lation  between  the  laminar  flame  speed  and  RH  to  hydrogen/ 
butane/air  mixtures.  Both  measurements  and  computations  show 
that  the  linear  correlation  still  holds  for  hydrogen/butane-air 
mixtures,  the  following  correlation  was  proposed: 

-  Rh)  —  Su((pj:-  0)1  k((pj:)Rij  (f3) 

As  we  mentioned  previously,  Eq.  (13)  simply  correlates 
the  laminar  flame  speed  with  the  parameters  (pF  and  RH,  while 
the  empirical  correlation  of  Eq.  (9)  is  more  complicated.  However, 
the  underline  factor  is  consistent:  the  laminar  flame  speed  of  the 
hydrogen/hydrocarbon  mixture  depends  strongly  on  the  mixture 
richness  and  hydrogen  addition  level. 

Unlike  the  statement  of  Yu  et  al.  [95]  that  the  slope  was  a 
constant  for  all  effective  equivalence  ratios,  Tang  et  al.  [103] 
discovered  that  the  slope  of  the  linear  correlation,  k,  depends 
strongly  on  the  effective  equivalence  ratio,  <pF,  as  shown  in  Fig.  24. 
The  linear  dependence  is  weakened  as  the  mixture  is  close  to 
stoichiometry.  In  other  words,  the  hydrogen  addition  effect  is 
minimized  at  around  stoichiometric  condition. 

Fundamentally,  when  a  certain  amount  of  hydrogen  is  added  to 
hydrocarbon,  the  laminar  flame  speed  is  influenced  by  three 
mechanisms:  (a)  the  thermal  effect,  caused  by  the  higher  adiabatic 
flame  temperature  of  hydrogen;  (b)  the  transport  effect  caused  by 
the  high  mobility  of  hydrogen;  (c)  the  kinetic  effect  due  to  the 
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Fig.  24.  Linearity  coefficient,  k,  quantifying  the  increase  in  flame  speed  with  RH 
[103]. 
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Fig.  25.  Sensitivity  coefficients  measuring  kinetic,  thermal,  and  diffusion  effects  on 
laminar  flame  speeds,  with  hydrogen  addition,  as  functions  of  RH  [103]. 

strong  reactivity  of  hydrogen.  Because  the  laminar  burning  flux 
(which  is  the  eigenvalue  for  flame  propagation)  is  given  by 

f2~Leexp(-Ea/R°Tad)  (14) 

in  which  the  Lewis  number  Le,  the  activation  energy  Ea,  and  the 
adiabatic  flame  temperature  Tad  can  be  considered  to,  respectively, 
represent  the  diffusion,  kinetic,  and  thermal  effects.  To  assess  the 
relative  importance  of  the  three  factors,  Eq.  (14)  was  differentiated 
to  give  the  overall  sensitivity  coefficient  as 

din  /  1  die  1  dTa  Ta  dTai 

d  In  Rh  2 Le  d  In  RH  2Tad  d  In  RH  2 T2ai  d  In  RH 

where  Ta=EaIR°  is  the  activation  temperature  and  the  three  terms 
on  the  right  respectively  represent  the  individual  sensitivity 
coefficients  for  diffusion,  kinetic,  and  thermal  effects. 

Fig.  25  shows  the  comparison  of  these  sensitivity  coefficients  as 
functions  of  RH  for  three  representative  lean,  stoichiometric  and 
rich  cases.  It  is  seen  that  for  all  tpF,  the  kinetic  term  is  the  largest 
sensitivity  factor  in  increasing  the  burning  flux.  The  sensitivity 
factor  for  the  thermal  effect  is  smaller,  but  nevertheless  still 
significant,  while  the  factor  for  the  diffusion  effect  is  substantially 
smaller  and  can  also  assume  negative  values  for  lean  mixtures. 
Furthermore,  all  three  sensitivity  factors  at  the  stoichiometric 
condition  assume  values  one  order  smaller  than  those  at  off- 
stoichiometric  conditions  of  tpF=0.&  and  1.4.  This  corroborates 
the  earlier  result  that  laminar  flame  speed  shows  the  weakest 
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dependence  on  RH  at  near  stoichiometric  condition  [103],  Further 
computations  of  other  hydrocarbons  showed  that  the  above  flame 
speed  dependence  on  RH  holds  for  methane,  ethane,  and  propane, 
indicating  that  there  exists  the  possibility  of  generalizing  the 
phenomena. 

5.2.  Flame  front  stability 

Accurate  determination  of  laminar  burning  velocity  is  always 
counteracted  by  the  presence  of  cells  over  the  flame  surface,  which 
increases  the  flame  surface  area  and  consequently  enhances  the 
flame  propagation  speed  [142],  especially  at  high  pressures  and/or 
for  lean  hydrogen  or  light  hydrocarbon  flames  and  rich  heavy 
hydrocarbon  flames.  The  flame  front  cellular  instability-induced 
self-acceleration  and  the  subsequent  self  turbulization,  and  even 
transition  to  detonation  [143],  could  potentially  promote  engine 
knock  [144,145].  For  this  reason,  cellular  instability  has  raised  a 
significant  amount  of  attention. 

Three  mechanisms  of  cellular  instability  are  well  established: 
hydrodynamic  instability  [146,147],  diffusional-thermal  instability 
[104,148,149],  and  instability  resulting  from  the  body  force  [150], 
The  body  force  induced  -  or  the  so-called  Rayleigh-Taylor 
instability  -  is  conspicuous  when  the  flame  speed  is  decreased, 
because  growth  rate  is  strongly  dependent  on  the  flame  speed 
[11  ].  Hydrodynamic  instability  was  first  recognized  theoretically  in 
1938  by  Darrieus  [146],  and  later  in  1944,  it  was  predicted 
independently  by  Landau  [147],  They  assumed  that  the  flame 
was  a  thin  interface  of  density  discontinuity  propagating  towards 
the  unburned  mixture  at  a  constant  speed,  the  growth  rate  of  the 
hydrodynamic  disturbance  is  proportional  to  the  thermal  expan¬ 
sion  ratio  and  thus  the  hydrodynamic  instability  is  strongly 
dependent  on  the  thermal  expansion  ratio  and  flame  thickness. 
It  should  be  noted  that  no  length  scale  was  embedded  in  the 
Darrieus-Landau  theory,  thus  the  flame  was  unconditionally 
unstable  to  all  wavelengths  of  hydrodynamic  disturbances.  Clanet 
and  Searby  [151]  established  a  perfectly  planar  laminar  flame  with 
an  acoustic  re-stabilization  technique,  as  shown  therein.  They 
were  among  the  first  to  experimentally  measure  the  growth  rate 
of  Darrieus-Landau  instability  and  their  observed  image  was  given 
in  Fig.  26. 

For  mixtures  in  which  the  deficient  reactant  is  also  the 
constituent  of  the  largest  diffusivity,  the  formation  of  cells  was 


Fig.  26.  Images  from  high  speed  film  of  growth  of  instability  [151 1.  Framing  rate  of 
500  images/s,  wave  length  of  2  cm  and  laminar  flame  speed  of  11.5  cm  /s. 


Fig.  27.  Images  of  fully  developed  cell  structure  for  rich  propane/air  flame  [148], 


observed  for  both  plane  flame  [148],  as  shown  in  Fig.  27,  and 
spherical  flame  [104],  as  shown  in  Fig.  28(b).  The  cellularity  in 
such  cases  is  a  consequence  of  a  diffusional-thermal  instability, 
which  results  from  the  competing  effects  of  heat  conduction  from 
the  flame  and  reactant  diffusion  towards  the  flame  [149,152,153], 
Lewis  number  (Le)  is  defined  as  the  ratio  of  heat  diffusivity  of  the 
mixture  to  mass  diffusivity  of  the  limiting  reactant.  When  Le  is 
below  some  critical  value,  Le*  (slightly  lower  than  unity), 
diffusional-thermal  instability  could  be  observed  during  the  initial 
phase  of  propagation,  i.e.,  the  flame  radius  is  the  order  of  flame 
thickness  [152], 

Obviously  for  mixtures  not  apart  from  stoichiometry,  the 
formation  of  cells  on  the  flame  surface  results  from  the  combined 
influence  of  the  hydrodynamic  and  diffusional-thermal  instabil¬ 
ities.  The  opposite  nonequidiffusive  behavior  of  heavy  hydrocar¬ 
bon  (propane,  for  instance)/air  and  hydrogen/air  was  well 
established  [138,154],  However,  there  are  limited  reports  on  the 
cell  formation  of  hydrogen  enriched  hydrocarbon  flames,  espe¬ 
cially  at  high  pressures  where  the  hydrodynamic  instability 
mechanism  is  conspicuous. 

Rayleigh-Taylor  instability  is  conspicuous  only  when  the  flame 
speed  is  very  low,  thus  the  effect  of  hydrogen  addition  tends  to 
prohibit  this  instability  mechanism  because  the  addition  of  hydro¬ 
gen  will  increase  the  laminar  flame  speed. 

Law  and  Kwon  [128]  studied  the  propane  substituted  hydro¬ 
gen/air  flame  at  atmospheric  pressure  and  found  that  the  flame 
front  instability  of  the  lean  hydrogen  air  flame  is  progressively 
moderated,  or  diminished,  as  SFHC  increases;  this  is  clearly  shown 
from  the  flame  images  (as  shown  in  Fig.  28(a)).  Other  results  of 
Law  et  al.  [135]  and  Tang  et  al.  [155]  were  also  shown  in  the  figure. 
The  propane  substitution  extent  parameter  SFHC  and  the  propane 
mixing  ratio  a  are  all  converted  to  the  simply  understanding 
hydrogen  percentage  XH2,  as  tabulated  at  the  bottom,  or  on  the 
right,  of  each  series  of  images.  It  is  expected  that  both  the 
diffusional-thermal  instability  and  the  hydrodynamic  instability 
were  prohibited  with  increasing  SFHC,  a  or  decreasing  XH2,  respec¬ 
tively  in  Ref.  [128,135,155],  due  to  suppression  of  diffusional- 
thermal  cells  and  the  substantial  increase  of  flame  thickness,  as 
more  propane  is  present  in  the  mixture.  The  Lewis  number  (which 
accounts  for  the  thermal-diffusion  effect)  is  conventionally  esti¬ 
mated  for  sufficiently  off-stoichiometric  mixtures  from  the  free 
stream  values  of  the  mixture  transport  properties.  For  blends  of 
fuel  with  different  transport  properties,  however,  the  Lewis 
number  should  be  a  global  flame  parameter,  evaluated  through 
flame  response  or  with  specific  weighted  averaging.  As  shown  in 
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Fig.  28.  Schlieren  photographs:  (a)  C3H8-substituted  H2-air  (0=0.6  at  NTP).  Ref.  [128].  (b)  Hydrogen,  hydrogen/propane  (propane  mixing  ratio  =0.25,  0.5  and  0.75)  and 
propane  in  air  at  5  atm.  Adapted  from  Ref.  [135].  (c)  Lean  (0=0.8),  stoichiometric  and  moderately  rich  (0=1.2)  propane/hydrogen  mixtures  at  5  atm  (<r,<5,  mm.Ka ).  Adapted 
from  Ref.  [155].  (d)  Adequately  rich  (0=1.4)  propane/hydrogen  mixtures  at  5  atm.  Adapted  from  Ref.  [155]. 


Fig.  29(a),  the  flame  Lewis  number  is  decreased  with  the  increase 
of  XH2  (decreasing  propane  mixing  ratio  a),  indicating  the 
thermal-diffusionally  de-stabilizing  effect  of  hydrogen  addition. 

For  the  rich  flames,  the  effect  of  propane  substitution  is  not 
clear  because  both  the  flame  thickness  and  the  density  ratio 
increase  substantially  as  the  SFHC  increases.  The  rich  hydrogen/ 
air  flames  are  thermal-diffusionally  stable  and  thus  Law  and  Kwon 
concluded  that  increasing  propane  would  not  likely  promote  flame 
front  instability;  however,  their  images  were  taken  at  atmospheric 
pressure,  where  the  hydrodynamic  instability  effect  is  not  as 
strong.  Another  series  of  Schlieren  images  by  Law  et  al.  [135] 
was  taken  at  higher  pressures,  but  no  results  of  fuel  rich  condi¬ 
tions  were  presented  therein. 

Tang  et  al.  [155]  showed  that  for  adequately  rich  propane/air 
flame,  the  net  effect  of  adding  hydrogen  to  flame  front  instability 
depends  on  the  competition  between  the  thermal-diffusion 


mechanism  and  the  hydrodynamic  mechanism,  and  they  did 
observe  that  the  flame  front  was  first  destabilized,  then  stabilized 
with  the  increase  of  the  hydrogen  addition,  as  shown  in  Fig.28(d). 
For  XH2  less  than  40%,  although  the  Lewis  number  increases  (as 
shown  in  Fig.  29(b)),  and  thus  the  thermal-diffusion  effect  tends  to 
stabilize  the  flame;  flame  front  instability  is  still  promoted  because 
of  the  hydrodynamic  destabilizing  effect  due  to  the  significant 
decrease  in  flame  thickness  (<5i).  However,  forXH2  larger  than  40%, 
the  Lewis  number  increases  more  significantly  and  the  flame  tends 
to  be  stabilized  as  XH2  increases.  This  implies  the  dominant 
influence  of  the  Lewis  number  compared  to  those  of  other 
parameters,  such  as  a,  <5|,  and  Ka.  The  non-monotonic  behavior 
at  fuel  rich  conditions  indicates  that  there  exists  a  mixture 
composition  (XH2=40%),  at  which  the  flame  is  the  most  unstable 
due  to  the  combination  of  effects  of  the  hydrodynamic  and 
diffusional-thermal  instabilities. 
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Fig.  29.  Global  flame  Lewis  number  for  propane/hydrogen/air  mixture  as  a  function 
of  XH2-  (a)  Evaluated  from  flame  response,  Refs.  [102,135].  (b)  Weighted  average 
value  [155]. 

6.  Effects  of  hydrogen  addition  on  lean  premixed  combustion 

Operating  an  engine  with  excess  air,  under  lean  conditions,  has 
significant  benefits  in  terms  of  increased  engine  efficiency  and 
reduced  emissions.  However,  under  high  dilution  levels,  a  lean 
limit  is  reached  in  which  combustion  becomes  unstable,  signifi¬ 
cantly  deteriorating  drivability  and  engine  efficiency  and  limiting 
the  full  potential  of  lean  combustion.  Due  to  the  high  laminar 
flame  speed  of  hydrogen,  adding  a  hydrogen-rich  mixture  with 
gasoline  into  the  engine  helps  to  stabilize  combustion,  extending 
the  lean  limit. 

Lean  premixed  combustion  (LPC)  is  also  conceptually  advanta¬ 
geous  in  gas  turbines  due  to  its  lower  flame  temperature,  which 
results  in  less  thermal  load  and  lower  emission  of  NOx.  However, 
LPC  is  susceptible  to  combustion  instabilities  because  of  the  weak 
burning  and  extinction  caused  by  local  stretching  and  heat  loss 
[156],  One  approach  to  solving  the  combustion  stability  problem  of 
LPC  is  through  the  addition  of  H2,  because  hydrogen  can  sustain 
vigorous  burning  at  relatively  lower  temperatures  compared  to 
hydrocarbon  [95], 

6.1.  Swirl  stablized  lean  premixed  combustion 

There  have  been  several  studies  of  the  effects  of  hydrogen 
addition  on  the  lean  stability  limit  of  swirl-stabilized  premixed 
flames  [157-159]  and  stagnation  flames  [30,156],  These  works 
demonstrated  a  significant  extension  of  the  lean  stability  limit  by 
the  addition  of  hydrogen. 

Strakey  et  al.  [159]  measured  the  lean  blowout  equivalence 
ratio  of  natural  gas  and  air  mixtures  with  a  hydrogen  addition  in  a 
swirl-stabilized  burner;  they  found  that  when  the  hydrogen 
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Fig.  30.  Lean  stability  limit  data  of  hydrogen  enriched  flames.  Data  obtained  from 
Strakey  et  al.  [159]  (swirl-stabilized  burner  for  H2+natural  gas,  squares  (□):  1  atm; 
circles  (°):  4  atm;  up  triangles  (*):  10  atm);  Schefer  [157]  (swirl-stabilized  burner 
for  H2+methane  at  1  atm,  down  triangles,  v),  and  Guo  et  al.  [30]  (counter-flow 
flames  of  H2+methane  at  1  atm,  diamond,  o). 

concentration  XH2  in  the  fuel  mixtures  is  increased  from  0  to 
80%,  the  lean  blowout  equivalence  ratio  is  reduced  from  (p~  0.46 
to  </>~0.30 .  Additionally,  their  experimental  data  indicated  that  at 
different  pressures  the  lean  stability  limit  is  insensitive  to  a 
variation  of  pressure.  Schefer  [157]  investigated  stability  charac¬ 
teristics  of  premixed,  swirl-stabilized  methane  air  flames  with 
emphasis  on  the  effect  of  a  hydrogen  addition.  It  is  seen  in  Fig.  30 
that  Schefer’s  measurement  is  slightly  lower  than  that  of  Strakey 
et  al.  This  might  be  caused  by  the  hydrocarbon  compositions  in  the 
natural  gas. 

Guo  et  al.  [30]  measured  the  extinction  limit  of  premixed 
counter-flowing  methane  air  flames  with  a  hydrogen  addition  of 
up  to  40%;  the  extinction  limit  was  also  included  in  Fig.  30.  It  is 
seen  that  Guo’s  measurement  is  lower  than  that  of  Schefer  and 
Strakey  et  al.,  which  is  reasonable  because  an  unstable  flame  could 
be  extinguished  by  decreasing  the  equivalence  ratio.  Overall,  the 
extension  of  the  lean  stability  limit  by  the  addition  of  hydrogen  is 
well  illustrated  by  Fig.  30.  The  extinction  strain  rate  of  hydrogen 
enriched  lean  premixed  methane  flames  in  highly  strained 
counter-flowing  fields  has  been  experimentally  and  numerically 
investigated  by  Jackson  et  al.  [160];  they  experimentally  observed 
that,  with  the  increase  of  hydrogen  concentration  in  the  fuel 
blends,  the  extinction  strain  rate  is  significantly  reduced. 


6.2.  Pollutant  emissions 

Hydrogen  addition  can  significantly  reduce  the  emission  of  CO. 
However,  one  concern  about  adding  hydrogen  into  hydrocarbon  is 
the  increase  of  flame  temperature,  which  may  result  in  higher 
emissions  of  NO*.  Schefer  et  al.  [158]  investigated  the  pollutant 
emissions  from  the  swirl-stabilized  lean  premixed  natural  gas 
combustion  with  and  without  the  addition  of  hydrogen.  Fig.  31 
shows  the  radial  profile  of  CO  and  NOx  emissions  at  the  axial  plane 
of  z=5.1  and  20.3  cm.  It  is  seen  that  in  the  central  region,  the  CO 
emissions  at  both  axial  planes  from  natural  gas  combustion  range 
from  15  to  100  ppm,  while  the  hydrogen  enriched  flames  produce 
CO  values  of  less  than  10  ppm.  Fig.  31(b)  shows  that  the  NOx 
concentration  from  hydrogen  enriched  natural  gas  combustion  is 
higher  than  that  from  pure  natural  gas.  This  is  because  the 
combustion  of  natural  gas  requires  longer  residence  distance 
than  the  hydrogen  enriched  flames  in  order  to  reach  steady 
NOx  concentration  because  of  faster  combustion  of  with  hydrogen. 
At  z=20.3,  steady  NOx  concentration  is  reached,  as  shown  in 
Fig.  31(d),  and  both  fuels  emit  similar  level  of  NOx.  These  results 
demonstrate  that  a  significant  reduction  in  CO  emissions  is 
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Fig.  31.  CO  and  NO*  concentration  profiles  in  the  radial  direction.  Measurements  at  plane  z=5.1  and  20.3  cm.  Adapted  from  Schefer  et  al.  [158]. 


realized  by  hydrogen  addition  without  adversely  affecting  the  NO* 
emissions. 


7.  Concluding  remarks 

Increasingly  unavailable  fossil  fuels  and  greenhouse  gas  emis¬ 
sions  from  transportation  vehicles  are  driving  investigations  of 
alternative  fuels.  Among  them,  hydrogen  offers  the  greatest 
potential  benefit  to  the  environment  and  energy  supply.  Hydrogen 
is  the  most  abundant  element  in  the  universe;  it  is  a  versatile 
energy  carrier  that  can  be  made  from  a  variety  of  primary  energy 
sources  such  as  natural  gas,  coal,  and  biomass,  and  non-carbon 
energy  sources  such  as  solar,  nuclear,  hydroelectric,  and  wind. 
Hydrogen  can  burn  with  high  efficiency  and  produce  essentially 
zero  emissions.  However,  attractive  as  it  is,  a  hydrogen  economy  is 
challenged  by  technical,  economical  and  infrastructural  barriers. 
For  instance,  there  are  significant  difficulties  associated  with 
hydrogen  storage  due  to  its  high  flammability  limits,  low  ignition 
energy  (which  causes  safety  problems),  and  its  low  volumetric 
energy  content,  which  requires  extra  energy  for  high  pressure 
storage.  Additionally,  there  exists  some  unresolved  issues  for  pure 
H2  combustion,  such  as  knock,  detonation,  pre-ignition,  and  flash¬ 
back.  For  these  reasons,  fossil  fuels  still  dominate  the  current 
primary  energy  supply.  To  bridge  this  situation,  a  more  rational 
approach  is  to  use  hydrogen/hydrocarbon  fuel  blends,  which 
synergistically  resolve  storage  and  combustion  problems  asso¬ 
ciated  with  burning  pure  hydrogen  or  pure  hydrocarbons. 

In  this  work,  recent  progress  in  the  fundamental  investigations 
on  the  combustion  characteristics  of  hydrogen/hydrocarbon  fuel 
blends  has  been  reviewed.  Experimental  results  show  that  the 
flammability  limits  of  hydrocarbon  are  extended  with  the  addition 
of  hydrogen.  Also  Le  Chatelier’s  Rule  can  predict  the  dependence  of 
a  lower  flammability  limit  with  reasonable  accuracy.  The  mini¬ 
mum  ignition  energy  of  hydrogen/hydrocarbon  blends  shows  a 
power-law  dependence  on  the  quenching  distance  and  both 
decrease  significantly  with  the  increased  hydrogen  fraction  XH2. 


Ignition  of  methane  and  air  stream  in  a  supersonic  mixing  layer  is 
promoted  with  hydrogen  addition  because  of  the  accelerated 
production  of  H  and  O  radicals.  Non-premixed  ignition  tempera¬ 
tures  in  counter-flowing  methane  vs.  hot  air  streams  are  decreased 
significantly  with  a  small  amount  of  hydrogen,  and  this  effect  is 
moderated  as  hydrogen  concentration  is  increased  to  7%  and  the 
effect  of  hydrogen  addition  effect  disappears  when  the  hydrogen 
fraction  is  higher  than  30%. 

Auto-ignition  of  hydrogen  enriched  hydrocarbons  by  shock- 
waves  has  been  reviewed  in  detail.  With  an  increase  of  the 
hydrogen  fraction,  ignition  behavior  of  hydrogen/hydrocarbon  fuel 
blends  in  Argon-diluted  oxygen  exhibits  three  distinct  behaviors. 
For  a  hydrogen  fraction  less  than  40%,  the  dependence  of  ignition 
delay  time  on  pressure  and  equivalence  ratio  resembles  that  of 
pure  hydrocarbon  and  ignition  is  dominated  by  methane  chem¬ 
istry.  For  a  hydrogen  fraction  higher  than  80%,  the  dependence  of 
the  ignition  delay  time  on  pressure  and  equivalence  ratio  resem¬ 
bles  that  of  pure  hydrogen  and  ignition  is  dominated  by  hydrogen 
chemistry.  For  a  hydrogen  fraction  in  the  range  between  40% 
and  80%,  ignition  delay  time  dependence  on  pressure  and  equiva¬ 
lence  ratio  resembles  neither  hydrogen  nor  hydrocarbon,  and 
ignition  is  governed  by  the  combined  chemistry  of  hydrogen  and 
hydrocarbon. 

The  laminar  flame  speed  of  hydrogen/hydrocarbon  fuel  blends 
is  found  to  exhibit  three  behavior  regimes.  When  XH2  increases, 
the  laminar  flame  speed  of  the  fuel  blend  increases,  but  the 
increment  is  small;  the  increment  becomes  significant  when  XH2 
is  higher  than  about  40%.  When  XH2  is  further  increased  to  more 
than  80%,  the  laminar  flame  speed  of  the  mixture  increases  even 
more  substantially.  When  a  specifically  defined  hydrogen  addition 
parameter  RH  is  defined,  the  flame  speed  of  n-butane/hydrogen  is 
found  to  increase  linearly  with  RH  due  to  altered  thermodynamics, 
diffusivity,  and  chemical  kinetics  with  hydrogen  addition.  Addi¬ 
tionally,  the  effect  of  hydrogen  is  minimized  at  around  the 
stoichiometric  condition.  Computation  with  other  fuels  such  as 
methane,  ethene,  and  propane  showed  that  the  laminar  flame 
speed  dependence  on  RH  is  generous  for  these  hydrocarbons. 
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There  are  three  mechanisms  that  influence  flame  front  instability 
behavior:  the  body-force  effect,  the  hydrodynamic  effect  and  the 
thermal-diffusion  effect. 

Hydrogen  addition  will  reduce  the  body-force  effect  because  of 
increased  flame  speed.  For  lean  heavy  hydrocarbons  (C3  or 
higher),  a  hydrogen  addition  will  decrease  the  overall  Lewis 
number,  resulting  in  a  promoted  thermal-diffusional  instability. 
While  for  sufficiently  rich  heavy  hydrocarbons,  thermal- 
diffusional  instability  tends  to  be  inhibited  by  the  addition  of 
hydrogen.  Hydrodynamic  instability  is  promoted  by  the  increase  of 
hydrogen  because  of  the  decreased  flame  thickness.  A  hydrogen 
addition  increases  the  resistance  to  the  strain  rate,  leading  to  an 
extended  lean  stability  limit  and  better  stabilized  swirl  burner 
combustion.  Additionally,  a  significant  reduction  in  CO  emissions 
is  realized  by  the  addition  of  hydrogen,  without  adversely  affecting 
the  NOx  emissions. 
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